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ABSTRACT 



Measurements of a 1230-km baseline were made during an eight-week period in the 
fall of 1987 using Trimble 4000SX single-frequency, five channel Global Positioning 
System (GPS) receivers. Twenty-eight days of carrier phase data were processed using 
correlated triple differences with fixed satellite orbits, the broadcast ephemerides, a 
modified Hopfield tropospheric model, and without ionospheric correction to determine 
the accuracies and precisions of the slope distance and baseline components. The data 
were processed in ever increasing observing sessions to determine the optimized obser- 
vation periods required to achieve various orders of geodetic accuracies. 

The accuracies of the slope distances were better than 1.0 ppm for any observing 
period. The day-to-day repeatabilities of the slope distance measurements were better 
than 1.0 ppm ( 2cr) after 20 minutes of observations. Accuracies and repeatabilities (2a) 
of the baseline components were better than 10.0 ppm after 20 minutes of observations. 
The correlated triple difference results were on the order of previous GPS surveys that 
used higher resolution differencing or external timing aids. Discussions include the ef- 
fects of ephemeris, tropospheric and ionospheric errors, and dilution of precision. 

Observation periods and mean slope distance errors were reduced when observations 
started close to and included the infinite peak of the Position Dilution of Precision 
(PDOP). The smallest variances were associated with observations about the infinite 
PDOP peak. 
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I. INTRODUCTION 



A priori knowledge of the observation periods required to achieve specified orders 
of geodetic accuracies is important in planning efficient and productive geodetic surveys. 
While terrestrial surveys have specified field and processing procedures and standards to 
categorize the geodetic accuracies of surveys [Federal Geodetic Control Committee, 
1984], only recently have standards been proposed for surveys conducted with the Global 
Positioning System (GPS) [Federal Geodetic Control Committee, 19S6]. Among the 
proposed requirements are standards for the length of observing periods and satellite 
geometry. 

Field studies by Remondi [1984] and numerical simulations by Fell [1980], Langley 
et al. [1984], and Landau and Eissfeller [1986] studied optimized observation periods, but 
for baselines less than 100 km. Cannon et al. [1985], Bock et al. [1984], Goad et al. 
[1985], Mader and Abell [1985], and Bertiger and Lichten [1987] conducted long baseline 
surveys, but did not study optimized observation periods. One of the objectives of this 
thesis is to fill the gap between the above studies, i.e., exannne the optimized observation 
periods for a long baseline. 

The optimized times will be examined using the correlated triple difference carrier 
beat phase observable because of its insensitivity to integer ambiguities and loss of lock 
of the GPS carrier by the receiver. Another objective of this thesis is to add to the body 
of triple difference accuracy testing following a recommendation by Remondi [1984, p. 
259]: “More testing is required to establish the full accuracy potential of the triple dif- 
ference method. ” 

GPS carrier phase and pseudorange measurements were made during an eight-week 
period in the fall of 1987 with Trimble 4000SX single-frequency, 5 channel receivers. 
The long baseline is approximately 1230 km in length between the National Oceanic and 
Atmospheric Administration's (NOAA) Western Regional Center located at Sand Point 
in Seattle, Washington, and the Naval Postgraduate School (NTS), Monterey, 
California. The baseline was determined by locating the positions of its ends by con- 
necting them by independent short baseline surveys from nearby Very Long Baseline 
Interferometry (VLB I) horizontal control points. The results of those surveys form the 
reference to which accuracy will be determined. 
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Additionally, studies for repeatability were conducted following another recommen- 
dation by Remondi [1984, p. 263] to enhance the capabilities of GPS measurements. 
The recommendation was to perform extensive repeatability studies on non-varying 
baselines for verifying and improving the GPS modelling. 
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II. BASELINE DETERMINATION USING GPS 



A. INTRODUCTION 

The Department of Defense's (DOD) Global Positioning System is intended to 
provide accurate positioning and precise timing for navigation purposes by broadcasting 
codes superimposed on two radio carrier frequencies from satellites. The satellites are 
placed in a constellation so that at least four satellites are visible globally. The Precise 
Code (P code) will be limited to authorized DOD users. The Coarse Acquisition (C A) 
code provides real-time accuracies to about 100 m [Baker, 1986] and is available to 
anyone. 

The codes provide their transmit times, satellite orbit and clock information, and 
information to enable any receiver to lock onto other GPS satellites. The orbital infor- 
mation (ephemeris) provides the position of the satellite. The receiver measures the time 
delay between the receipt of the C/A code and its transmission time. The time delay can 
be transformed into an apparent slant range from the satellite's known position to de- 
termine the location of the receiver. Since it includes delays due to receiver clock errors 
and the effects of atmospheric refraction, the apparent slant range is referred to as the 
pseudorange. A minimum of four satellites are required to solve the system of range 
equations for the receiver's coordinates and clock errors. 

While the C'A code provides real-time location, it does not meet the accuracy re- 
quired for precise geodetic work. Nevertheless, GPS makes possible a higher resolution 
via the carrier signal. Though the carrier itself does not contain the orbital and timing 
information, which would have to be supplied by some other means, it does offer a 
higher resolution because of its 19-cm wavelength. 

B. THE MONTEREY-SAND POINT BASELINE 
1. General 

The length of the Monterey-Sand Point baseline (Figure 1) was computed by 
differencing the World Geodetic System 1984 (WGS84) [Defense Mapping Agency, 
1987] Cartesian coordinates determined for Monterey and Sand Point by short baseline 
GPS surveys from known horizontal control points. The precision and agreement with 
terrestrial survey results of short baseline GPS measurements using single frequency, 
double difference solutions are well documented ( e.g ., [Remondi, 1984], [Goad and 
Remondi, 1984], [Bock et al., 1984]). 
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Figure 1. Monterey-Sand Point Baseline and environs: Insets not drawn to scale. 

At Sand Point, on-site meteorological measurements were made near the middle 
of the observing session. For the Monterey antenna determination meteorological 
measurements were made every half hour and the mean of all the measurements was 
used in the processing. The Trim640 solutions were obtained using uncorrelated double 
differences and estimating initial integer ambiguities. The ambiguities were fixed to the 
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integer values that produced the smallest residuals. A tropospheric factor was estimated 
along with the integer ambiguities and the baseline components in the least-squares 
processing. 

The horizontal control points used for the reference stations in determining the 
coordinates of the antennas were mobile Very Long Baseline Interferometry (VLBI) 
sites. The NAD83 Cartesian coordinates for the VLBI sites were provided by the the 
Gravity, Astronomy and Space Geodesy Branch of the National Geodetic Survey (NGS) 
[Abell. 1987]. The XAD83 coordinates were determined in August 1987 from a global 
adjustment of VLBI surveys. Carter et al. 11985] described the NGS VLBI program. 
The Defense Mapping Agency Hydrographic'Topographic Office (DMAHTC) validated 
the direct transformation of the VLBI Cartesian coordinates to WGS84 Cartesian coor- 
dinates [Kumar, 1988]. 

2. Monterey Coordinates 

The coordinates for the Monterey antenna location were determined by aver- 
aging two surveys conducted on separate days from the VLBI site FT ORD XCMN 
1981. Table 1 lists the observing sessions used to determine the WGSS4 coordinates of 
the Monterey antenna. Table 2 lists the results of the GPS surveys. 



Table 1. MONTEREY ANTENNA LOCATION SURVEYS: From FT ORD 

NCMN 1981. PST = Pacific Standard Time. 



Date 


Start 

Time 

(PST) 


End 

Time 

(PST) 


Number of 
Double 
Difference 
Observations 


Percent 

Rejected 


Slope 

Distance (nr) 


a b 


09 16 87 


1018 


1210 


996 


7 


12139.S125 


0.025 


09 18 87 


0914 


1120 


1136 


4 


12139.8119 


0.025 



The Monterey coordinates were computed from: 

Cmry ~ Cord + \ 

where: 

C MRr Monterey coordinate 

C 0RD Ft Ord coordinate 

A c Variance-weighted mean baseline component 

The uncertainties of the Monterey coordinates, a c , were computed from: 
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Table 2. RESULTS OF MONTEREY ANTENNA LOCATION 
SURVEYS: WGS84 Cartesian Coordinates (meters). 



DATE 


FT ORD X 


c ORD X 


AX 


°ZX 


MRY X 


G MRY X 


09/16 87 
09 18.87 


-2697026.493 

-2697026.493 


0.007 

0.007 


-10313.601 

-10313.602 


0.032 

0.037 


-2707340.093 

-2707340.095 


0.032 

0.037 


Mean: -2707340.094 


0.036 


DATE 


FT ORD Y 


° ORD Y 


AY 


°*Y 


MRY Y 


G MRY Y 


09 16,87 
09' 18 87 


-4354393.309 

-4354393.309 


0.010 

0.010 


917.693 

917.689 


0.048 

0.049 


-4353475.617 

-4353475.620 


0.049 

0.050 


Mean: -4353475.618 


0.051 


DATE 


FT ORD Z 


°ORD Z 


AZ 


G *z 


MRYZ 


° MRYZ 


09/16.87 
09 18/87 


3788077.778 

3788077.778 


0.009 

0.009 


-6337.391 

-6337.389 


0.041 

0.044 


3781740.387 
37S 1740.389 


0.042 

0.045 


Mean: 3781740.388 


0.044 



, / 2 .2 

Or = + \ Or OT 

C MRY — \ c ORD \ 

where a c is the uncertainty of the Ft Ord coordinate. 

One month prior to the surveys originating from Ft Ord, two other surveys were 
conducted from the satellite Doppler horizontal control point NAVAL POST GRAD 
3 1965,'DOPPLER. The Doppler station is approximately 300 m north of the Monterey 
antenna location (Figure 1). The double differenced GPS carrier phase solutions of the 
two Doppler-originating surveys yielded the following Monterey coordinates: 

X -2707339.725 m 

Y -4353475.654 m 

Z 3781740.264 m 

The baseline components had an uncertainty of ±0.002 m, but the Doppler station has 
an uncertainty of ±2 m in each coordinate before transformation to WGS84. The results 
of the Doppler surveys were not used in determining the Monterey coordinates because 
of the large uncertainty in the Doppler station location. The three-dimensional positions 
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ol the Monterey antenna from the Doppler and the VLBI originating surveys agree to 
better than 0.5 m. 

The two days of pseudorange observations at the Monterey receiver were each 
subjected to a least-squares estimation and then averaged together to yield the Monterey 
coordinates: 

X -2707334.248 ±1.3 m 

Y -4353466. S07 +0.6 m 

Z 3781741.330 ±2.2 m 

where the two-day standard deviations are given. The deviation of the pseudorange 
from the differenced carrier phase derived Monterey coordinates is expected because of 
the coarser resolution of the C/A code and because the pseudoranges are corrected for 
neither tropospheric nor ionospheric delays. 

3. Sand Point Coordinates 

The Sand Point antenna coordinates were determined by one 90-minute GPS 
survey (Table 3) from the mobile VLBI site Aviation 2 which is 530 m distant 
(Figure 1). 



Table 3. RESULTS OF SAND POINT ANTENNA LOCATION 



SURVEY: From double difference carrier phase solutions. Distances 

and WGS84 Cartesian Coordinates are in meters. 



Aviation 2 Coordinate 


g a n 


A 




Sand Point Coordinate 


a SFr 


X= -2295347.760 


0.017 


-408.361 


0.002 


-2295756.121 


0.017 


Y= -363S029.429 


0.028 


330.200 


0.002 


-3637699. 22S 


0.02S 


Z= 4693408.964 


0.032 


73.813 


0.003 


4693482.777 


0.032 



The Sand Point antenna coordinates computed from the pseudorange data 



were: 

X -2295749.623 m 

Y -3637694.279 m 

Z 4693488.472 m 

and standard deviations are not given because this is a single observing session and the 
program does not provide solution standard deviations. 
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4. Monterey-Sand Point Baseline Components 

First the Monterey-Sand Point baseline components were computed by sub- 
tracting the Sand Point Cartesian coordinate from the respective Monterey Cartesian 
coordinate. The uncertainties in the baseline components were computed as the square 
root of the sum of the squares of the uncertainty of the Monterey and Sand Point co- 
ordinates. The slope distance was computed as the square root of the sum of the squares 
of the baseline components. The slope distance uncertainty was computed as: 



V ( 



' AX 



A) 



AZ 



'AX 



x) + ( u a A>) + ( l a Az) 



where b is the slope distance. Using the information from Tables 2 and 3 gives 
Table 4, which is used as a standard to estimate accuracy. 



Table 4. STANDARD BASELINE DISTANCES AND 2-SIGMA VALUES. 



AX 


-411583.973 


± 0.080 m 


(0.19 ppm) 


AY 


-715776.390 


± 0.1 16 m 


(0. 16 ppm) 


AZ 


-911742.388 


± 0.109 m 


(0.12 ppm) 


Slope Distance 


1230045.280 


± 0.109 m 


(0.09 ppm) 



C. THE ONE WAY CARRIER BEAT PHASE 

The development of the carrier beat phase technique and a model for its application 
are given by Remondi [1984]. The phase measurement is done by beating the received 
carrier with a local oscillator internal to the GPS receiver. The slant range from a GPS 
receiver to a GPS satellite can be modelled in terms of the time it takes the signal to 
travel or the number of cycles that occur between the satellite and the receiver. The 
range in cycles will consist of an integer and fractional number of cycles. When a GPS 
receiver locks onto the carrier signal, it can immediately measure the fractional part and 
begin counting subsequent integer cycles, but it cannot measure or account for the initial 
integer number of cycles that preceded the initial fractional part. These missing cycles 
bias subsequent measurements and are called the initial integer ambiguity biases. 

The signal does not take a direct path to the receiver as it is refracted by the 
ionosphere and troposphere. Additional errors are caused by the satellite deviating from 
its predicted orbit, errors in the satellite clock, and errors in the receiver clock. Anti- 
cipating that the observables will be used in relative positioning, that a single frequency 
receiver will be used and ignoring other error sources, such as multipath, the one-way 
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carrier beat phase, <t> b (r,s.i), observed at receiver, r, from satellite, j, at observation 
epoch, /, can be modelled to first-order as: 



where: 



f 

<t> b {r,s,i) = 0X0 - <t> r {r,i) - -j- p(r,s,i) 

+ CMl Ifs +fr ~ 4 0('',V')] 
- t(r ,/) + /l(r,5,l) 

/ Epoch identifier 

r Receiver identifier 

s Satellite identifier 

0X0 Phase of received carrier signal 

Phase of receiver generated carrier signal 
f Transmitted frequency of carrier signal 

p{r,s,i ) Satellite-receiver slant range 

p(r, 5 .i) Time rate of change of slant range 

c Speed of light 

f Receiver generated carrier frequency 

r (r,i) Tropospheric delay 

.-f 1 ) Initial integer ambiguity 

ai,o=^(/)-44 

C ( 2,0 = {(0 + 



<5(0 



c(i) Mean clock offset for both receivers 

(5(0 Clock drift between both receivers 



( 2 . 1 ) 



and parentheses do not indicate factors or functions, but simply enclose identifiers. 
Brackets indicate factors. 



D. DIFFERENCING THE ONE-WAY CARRIER BEAT PHASE 
1. Single Difference 

The single difference ( SD ) is formed by differencing carrier beat phase observa- 
bles from two receivers at the same observation epochs. Following Remondi [1984], 
taking the difference, expanding the £(r,0 terms of Equation (2.1), and expressing dif- 
ferences in the between-satellite and between-receiver phases as f s d(i) gives: 
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SD{s,i) = <j> b (2,s.i ) - <J> b ( 1 ,s,i) 

= 4Cp(2.w)-p(1a«*)] 

+ 4 Cp( 2.5,/) - p( l, J, /)]<!(/) ^ 

— C/t — /2^<?( ; ) + /s<H') 

+ 4“ Cp( 2,5,/) + p( 1,5, /)]<)(/) 

+y;[r(2,/) - + 7(2,s,l) - 7(l,s,l) 

where the terms are the same as Equation (2.1). Single differencing reduces or eliminates 
satellite orbital and clock errors because they are common to both receivers. 

2. Double Difference 

A double difference (DD) is formed by differencing single differences between a 
reference satellite, h, and another satellite at the same epoch: 

DD(h,s,i) = SD{h,i) - SD(s,i) 

Because the differences at the same epoch are taken with the same reference 
satellite, the double differences for each epoch are correlated. The advantage of the 
double differences is that the clock dependent terms - f s 5(i) and [f 2 - are elim- 

inated. The significance of the removal of those terms is to reduce from nanoseconds 
to microseconds the timing accuracy required to achieve one cycle accuracy. The 
Trimble 4000SX achieves sub-microsecond accuracy by using the C/A code timing in- 
formation [Ashjaee, 1 985]. 

3. Triple Difference 

A triple difference (77)) is formed by differencing the double differences for the 
same satellite pair at some integer number of succeeding epochs, / : 

TD(i) = DD(h,s,i + 0 - DD(h, s,i) 

The advantage of the triple difference is that it eliminates all the time inde- 
pendent terms, namely the initial integer ambiguities, A(r,s, 1), and becomes insensitive 
to the initial ambiguities and any cycle slips when the receiver loses lock. 

The disadvantages of the triple difference are: another level of correlation, loss 
of resolution, and reduced number of observations. The triple differences are already 
correlated with respect to satellite because of the underlying double differences, and are 
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further correlated with respect to time because consecutive triple dilferenced observa- 
tions will have the DD{h,s,i + f) term in common. 

For short baselines, integer ambiguities can easily be resolved because unmod- 
elled errors are highly correlated between the two antenna sites and are mostly elimi- 
nated by the differencing. Algorithms can take advantage of the integer nature of the 
initial ambiguities and solve for them. At longer baselines the unmodelled errors are not 
as highly correlated and not eliminated by the differencing. These errors fold into the 
initial ambiguities, so that the ambiguities are no longer integers. In some cases the 
ambiguities cannot be resolved (e.g., Henson and Collier, 19S6, Tables 1 and 2). 

Because of the advantages the triple difference offers for long baselines, 1 use the 
triple difference scheme. While the triple difference can be decorrelated by forming a 
weight matrix [Remondi, 1984], only the correlated triple difference software was avail- 
able to me. 

E. ERROR EFFECTS 

For long baseline GPS surveys, the primary errors are satellite orbit errors, 
ionospheric and tropospheric delays [Remondi, 1984, and Beutler et al . , 19S6]. Orbit 
(ephemeris) errors are the result of the departure of the satellite from the broadcast 
ephemeris orbit. The ephemeris is a predicted orbit for the satellite. Orbit errors prop- 
agate directly into the baseline measurements when the GPS orbit coordinates are fixed 
in the processing [Hothem and Williams, 1985]. Orbit errors can be the dominant error 
source affecting the repeatability of long baseline measurements [Lichten and Border, 
1987], 

The magnitudes of baseline errors increase with increasing baseline length because 
of the increasing projection of the ephemeris error onto the baseline component [Fell, 
19S0]. Estimates of the broadcast ephemeris error range from 25 m [Beutler et al., 1986] 
to 100 m [Wells et al., 1986], The magnitude of the effect of the ephemeris error on 
baseline accuracy has been traditionally approximated [Beutler et al., 1984, Equation 
( 2 . 1 )]: 




where: 

b baseline length 

p slant range (receiver to satellite) 
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E b error in baseline length 

£ , error in slant range 

The slant range to a GPS satellite is about 20,000 km, which translates to a baseline er- 
ror ranging from 1 ppm to 4 ppm using Equation (2.3). 

It is expected that the ephemeris error for the Monterey-Sand Point baseline will be 
towards the lower end of the range because the ephemerides are uploaded prior to the 
satellites entering the Yuma Proving Ground (Russell and Schaibly, 1980] near the 
California border. The ephemeris linearization error specification is to 1 m per day 
(Wells et ai, 1986]. 

The ionosphere disperses the code (the group velocity) from the carrier phase (phase 
velocity) because the C/A code has a frequency of 1 MHz while the carrier signal upon 
which it is superimposed has a frequency of 1575 MHz. The effect is to increase the 
pseudoranges, but decrease the carrier phase derived ranges (Smith, 1987]. Field exper- 
iments (Beutler et al., 1986] showed that ionospheric dispersion shortens baselines on the 
order of a few tenths to perhaps 2 ppm. 

Ionospheric error is proportional to the Total Electron Content (TEC) along the 
signal path and the cosecant of the elevation angle of the satellite [Smith, 1987]. Thus 
the error is greatest for low elevation angles and least at the zenith. Wells et ai. (1986] 
estimated that the range errors due to the ionosphere are from 150 m at the horizon to 
50 m at the zenith. 

Ionospheric activity is a function of latitude, longitude, time of day, season, and 
sunspot activity. Ionospheric activity increases towards the equator and towards the 
sunlit portions of the earth. Diurnally, it has a minimum near 0600 local time with a 
maximum around 1600 local. Ionospheric activity increases with the peak in the sunspot 
cycle. The minimum in the current 1 1 -year sunspot cycle occurred in 1986. Upon these 
systematic characteristics sporadic ionospheric disturbances are superimposed. [Henson 
and Collier, 1986] 

The tropospheric error is proportional to the refractivity along the satellite-receiver 
path and proportional to the cosecant of the elevation angle [Martin, 19S0]. The index 
of refraction, n, is the ratio of the speed of light in a vacuum to the speed of light in a 
particular medium, in this case the troposphere. Because the refractive index is a small 
fraction greater than 1.0, a more convenient unit to work with is refractivity, A r , where 
A' = {n — 1) x 10 6 . The magnitude of the tropospheric biases range from 20 m for 10° 
elevation angles to 2.3 m at the zenith (Wells et al., 1986]. 
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Henson and Collier [1986| have shown that triple difference measurements are un- 
affected by path-dependent ionospheric bias errors, but by path-dependent gradient 
ionospheric errors, and Martin (19S0| has estimated that the combined ionospheric and 
tropospheric gradient errors are on the order of meters per hour and are proportional 
to the cosecant of the elevation angle. 
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III. DATA COLLECTION, PROCESSING, AND ANALYSIS 



A. EQUIPMENT CONFIGURATION 

1. Hardware 

A complete description of the Trimble 4000SX receiver is given by Trimble 
Navigation [1937a]. NPS operates three Trimble 40U0SX GPS Surveyor receivers, of 
which two were used in this study. The 4000SX is capable of observing the C A code, 
integrated Doppler and carrier beat phases of up to five satellites simultaneously. Its 
ability to use the C/A code allows the receiver to be used as a stand alone navigation 
system which determines position using Doppler-smoothed pseudoranges and velocities 
[Ashjaee, 1985], 

For precise relative positioning the 4000SX can transmit its data through an 
RS232 port to a microcomputer for storage on floppy disk for post-processing. The 
4000SX's ability to use the C/A code allows it to decode the GPS navigation messages 
so that it can track satellites automatically and determine. Most importantly, it uses the 
C/A code in a time transfer mode to determine any offset and drift of its own clock and 
thus provide accurate time tags for the observations without an external atomic clock 
or synchronization with the receiver at the other end of the baseline. 

The receivers were left on continuously to allow unattended data collection. 

Multipath-resistant Trimble microstrip antennas were installed at both the 
Monterey and Sand Point locations. 

2. Software 

For relative positioning, the receiver is controlled from the microcomputer by 
Version D of Trimble's Datalogger program. The reference position (the geodetic coor- 
dinates of the antenna) and the particular options chosen must be entered into the re- 
ceiver via the receiver keypad. The receivers were set to use the reference position height 
for point positioning when less than four satellites were available. 

Each observation session was initialized to log data when a minimum of four 
satellites were 15° above the antenna's horizon. Five satellites were designated for each 
observing session. The software logs the observables and receiver clock parameters to 
a floppy disk every 15 seconds and the C/A code-determined antenna position and Po- 
sition Dilution of Precision (PDOP) every five minutes. The GPS navigation message 
is logged to a separate file at the beginning of the session. 
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B. SATELLITE OBSERVATION PLAN 
1. Satellite Selection 

The same five satellite (SV) (6, 9, 11, 12, and 13) were used lor the entire eight- 
week observation period. These five satellites were visible at both stations for over 100 
minutes, and of these five satellites four were visible for three hours (Figure 2). 



START DATE/TINE: 1987/10/17 15:21:15. DAY OF YEAR 290 TOH 575675. 

STOP DATE/TIHE: 1987/10/17 18:16: 1. 0AY OF YEAR 290 TOH 584160. 

DATA AVAILABLE: 

STATION: SAND POINT 



SV 6 1 | 

SV 9 | | 

SV 11 | | 

SV 12 | | 

SV 13 I | 

STATION: MONTEREY 

SV 6 | | 

SV 9 I I 

SV 11 | | 

SV 12 | | 

SV 13 | | 



Figure 2. Satellite availability: Each dot represents 10 observations; each col- 

umn, 10 epochs. From Trim640 output. 



2. Position Dilution of Precision (PDOP) 

Trimble Navigation [1987b] recommends that observations include the time that 
the PDOP goes to infinity. The Federal Geodetic Control Committee [1986] notes that 
initial results from investigations indicate that best results may be achieved when the 
Geometric Dilution Of Precision (GDOP) is changing value during the observing ses- 
sion, and proposes that observing sessions start with a high GDOP and stop with a low 
GDOP. 

Position Dilution of Precision is a component of the GDOP. GDOP is a 
measure of how satellite geometry' degrades point position accuracy [Jorgensen, 1984]. 
For computational ease in the navigation solution the GDOP is defined as the square 
root of the trace of the covariance matrix of the errors in position and time with the range 
errors set to one [Milliken and Zoller, 1980]. The role and definition of GDOP and 
PDOP in GPS point positioning were applied to GPS relative positioning, i.e., good 
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PDOP would provide better accuracy than poor PDOP [King et al. 1985). Good and 
poor PDOP are shown in Figure 3. 

Landau and Eissfeller [1986], using numerical simulations in which they assumed 
a full 18 satellite constellation and a receiver that could track those satellites that mini- 
mized GDOP, found that better accuracy for triple difference solution for a 68 km 
baseline generally corresponded to high GDOP. They used a more complete GDOP that 
included consideration of ionospheric, tropospheric and satellite position errors which 
are neglected in the conventional GDOP. 

The 4000SX does not record GDOP, but it does record PDOP every live min- 
utes. PDOP relates to GDOP as: GDOP 2 = PDOP 2 4- TDOP 2 , where TDOP is the Time 
Dilution Of Precision, the error in the user clock bias multiplied by the velocity of light. 
The expected uncertainty in a GPS point positioning solution is a product of the PDOP 
and the expected slant range error. The difference between PDOPs at Monterey and 
Sand Point remained less than 1.0 for the entire eight week observing period. The PDOP 
peaks at two times in an observing session (Figure 4) - 60 minutes and 150 minutes. 
The PDOP peaks occur near when the satellites lie in a common plane causing the sol- 
ution of the linearized range equations to diverge and the PDOP becomes infinite 
[Jorgensen, 1984). 

C. METEOROLOGICAL PARAMETERS 

Meteorological parameters are needed for the tropospheric correction model used in 
the processing software. On-site meteorological observations were not available for 
Sand Point. Instead observations from the Weather Service Office (WSO) of the 
Seattle-Tacoma International Airport (SEA-TAC) were used for Sand Point. SEA-TAC 
is approximately 27 km from Sand Point (Figure 1). The NPS Department of Meteor- 
ology routinely collects real-time hourly observations of sea-level pressure, temperature 
and relative humidity from the National Weather Service's data network. Observations 
that pertained to GPS observing sessions were entered into a file that is accessed by the 
batch file building program. 

While the Federal Geodetic Control Committee [1986] has proposed the use of on- 
site meteorological parameters, researchers have had success using standard atmosphere 
parameters for satellite geodesy [Fell, 1976] or extrapolated meteorological data 
[Rathacher et al., 1986]. 
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Figure 3. Poor anti good PDOP: From King et al. (1985, Fig. 3.2). 




Figure 4. PDOP versus time 
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Hourly observations for Monterey were obtained from the NTS Department of 
Meteorology. The instruments were located approximately 300 m south of the 
Monterey antenna location (Figure 1 ). 

D. PROCESSING SOFTWARE 

A complete description of the Trimble supplied Trimvec software can be found in 
Trimble Navigation [1987b]. The data was processed using the Trimvec Trim640 pro- 
gram, Revision AB. Trirn640 limits processing to 700 epochs, so the first 700 epochs for 
each observing session are used. Sand Point was used as the reference station and its 
coordinates kept fixed in the least-squares processing. Sand Point was chosen as the 
reference station because four satellite availability occurred later at Sand Point than at 
Monterey. This avoided having to load the not-in-common epochs from Monterey at 
each processing. 

Trim640 uses the C/A code derived positions obtained at the lowest PDOP for the 
initial estimates of the baseline components. Trim640 culls the best C/A code position 
during the data loading. No ionospheric correction is provided by the software, and only 
the broadcast ephemeris can be used to compute fixed orbit satellite positions. In the 
triple difference processing, the only parameters estimated by the least-squares process- 
ing are the baseline components, AX, AY and AZ. 

A modified Hopfield tropospheric model [Goad and Goodman, 1974] is used to 
correct the carrier phase delay caused by the troposphere. The correction is a function 
of the atmospheric refractivity computed from surface meteorological values of pressure, 
temperature, and humidity, and the elevation angle of the satellite. Larger corrections 
are required for low elevation angles, as the signal travels a longer path through the 
troposphere. The model corrects for at least 90% of the tropospheric delay [Remondi, 
198 I]. Trim640 allows only one pressure, temperature and humidity entry for each site 
per session. 

E. PROCESSING PROCEDURES 

1. General 

To study optimized times, the data from each observing session were segmented. 
Each successive segment contained 10 minutes more data than its predecessor. For ex- 
ample, for an entire observing session that started when four satellites were available and 
stopped when less than four were available provides 175 minutes of observations. The 
first segment will use the first 10 minutes of data , the sixth segment will use the first 60 
minutes, while the eighteenth will use all 175 minutes. For each segment, the entire 
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processing was restarted from the data loading. Reloading the data for each segment 
takes considerably longer than using the Trim640 option to flag data for processing, but 
reloading was done so that the processing does not use a best C/A code position from 
later in the observing session. 

Convergence of the least-squares solution was achieved by doing five iterations 
using every tenth triple dillerence formed from every tenth double difference, followed 
by five iterations decreasing the triple and double difference increments to five, and 
finally five iterations using all triple differences formed from all double differences. 

Trim640 rejects those observations whose residuals exceed a multiple of the 
mean residual. The multiple of the mean residual is known as the edit multiplier. 
Trim640 uses 3.5 as the default value for the edit multiplier. I used the default value for 
the initial processing. 

Any segment that had more than ten percent of its observations rejected or 
whose solution slope distance standard deviation (a,) was greater than 10 m was re- 
processed. The reprocessing was identical to the initial processing except that before 
invoking the triple difference process the pseudoranges for both stations were subjected 
to separate least-squares adjustment. The pseudorange processing improves the C/A 
code derived initial estimates for the baseline components and corrects the carrier beat 
phase time tags. The carrier beat phase time tags are computed from the C/A code 
times, and are earlier than the C A code times. 

If the pseudorange processing failed to lower the rejections to ten percent, the 
edit multiplier was increased until the rejections reached ten percent. A ten percent re- 
jection level was observed for a few sessions and always occurred within the first thirty 
minutes of observations. The data were transferred to the Naval Postgraduate School's 
IBM 3033 computer for analysis. The data were analysed and graphics produced using 
the A PL-based GRAFSTAT program. 

To study the effects of reducing observation time, five case studies were under- 
taken in which the observation start times were changed for processing. Each case study 
followed the processing procedures outlined above. 

2. Batch Processing 

Processing is performed in a batch mode. A batch file passes parameters to a 
template. Trimble supplies command files that tell the Trim640 to use the template pa- 
rameters in processing. 
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The batch file is built using the program Batbld (Appendix A). Baibld builds a 
batch file by providing the appropriate file names and start and stop times. Baibld 
computes the appropriate meteorological parameters for each segment by locating the 
applicable weather observations from the weather observation file, interpolating values 
at the start time, computing running means from each hourly weather observation, then 
interpolating the running means to the stop time of each segment. Two millibars were 
subtracted from the SEA-TAC sea-level pressure to compensate for the 20-m elevation 
above sea-level for the Sand Point antenna. 

Initially, processing was done on an IBM XT with a math coprocessor and a 
hard disk. Processing the 18 segments of an observing session took ten hours of com- 
puting time. Later, processing was performed on an 80286 based microcomputer run- 
ning at 10 MHz with an 80287-8 math coprocessor that reduced the processing time to 
three hours. 

Two minor problems with Trim640 were discovered during the processing. 
First, large values in range differences were found when using the range differences 
rather than the pseudoranges to improve the C/A code positions. The data were for- 
warded to Trimble Navigation for evaluation and an error was found in their software. 
The error had no apparent effect upon carrier phase difference processing. Second, 
Trim640 is incompatible with one or more of the AST Research, Inc. device drivers 
supplied with the 80286 microcomputer: ram disk, print spooler, and extended memory. 
Removing the drivers allowed Trim640 to execute normally. 

At the conclusion of the batch processing, the slope distance, the baseline 
components, their standard deviations (ct s ), the number of observations, the number of 
observations rejected, and the RMS cycle fits were extracted from the Trim640 output 
file and collected into files that held the data for a particular segment for each case study. 

F. ANALYSIS PARAMETERS 

The statistical parameters that will be used to evaluate the results are the error, the 
sample mean and the standard deviation defined as [Davis et al ., 1981]: 

E d =C d —C s (5.1) 
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where: 

Ej Error for day d 
C d Measured component for day d 
C s Expected values 
E Mean error or bias 
a 2 Variance 

a Standard deviation of E d 

c;j Standard deviation of E 

Accuracy describes the closeness between the measurements and the expected values 
[Davis et al., 1981]. The degree of accuracy is determined to the magnitude of the mean 
error (£). The repeatability of the measurements will be expressed in terms of 2a because 
it approximates the 95% confidence level for single-dimension measurements [Federal 
Geodetic Control Committee, 1986]. The slope distance and individual baseline com- 
ponents are one-dimensional measurements. 

G. DATA AVAILABILITY 

Observations were made simultaneously at Monterey and Sand Point for an eight- 
week period beginning 29 September 1987. Observations were made Tuesday through 
Saturday except the days after federal holidays. Forty observing sessions were con- 
ducted of which 2S were used in the analysis and are listed in Table 5. The remaining 
12 days of observations were not used in the analysis for various reasons, which are 
listed in Table 6. 

For brevity the observing days will be referred to by their Julian day. Times in 
Table 5 are given in Pacific Standard Time (PST) rather than Universal Coordinated 
Time (UTC) for ease in the later discussions on diurnal effects. 
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Table 5. DAYS USED IN THE DATA ANALYSIS 



Date 


Julian Day 


Start Time 
(PST) 


End Time 

(PST) 


Number of 
Triple Difference 
Observations 


09 29 87 


272 


0833 


1129 


1599 


09 30 S7 


273 


OS 30 


1125 


1485 


10 01 87 


274 


0825 


1121 


1501 


10 02, 87 


275 


0822 


1118 


1479 


10 03 87 


276 


0817 


1113 


1518 


10 06.87 


279 


0805 


1100 


1563 


10 07 87 


2S0 


0801 


1056 


1490 


10.0S.S7 


2S1 


0757 


1053 


1508 


10 09 '87 


282 


0753 


1049 


1500 


10 10 87 


283 


0749 


1044 


1522 


10 14 87 


287 


0732 


1028 


1524 


10 17 87 


290 


0721 


1016 


1504 


10/20/87 


293 


0707 


1002 


1540 


10 21 S7 


294 


0703 


0958 


1562 


10/22/S7 


295 


0658 


0953 


1551 


10/24.87 


297 


0651 


0946 


1493 


10 27 87 


300 


0638 


0934 


1567 


10 28/87 


301 


0634 


0930 


1516 


10/29/87 


302 


0631 


0926 


1681 


11,03/87 


307 


0613 


0908 


1848 


11 04 87 


308 


0605 


0901 


1553 


11/05,87 


309 


0601 


0S57 


1695 


1 1/ 10 87 


314 


0541 


0837 


1535 


11/11/87 


315 


0536 


0832 


1927 


11/13/87 


317 


0528 


0824 


1923 


11/14 87 


318 


0524 


0819 


1550 


11,21/87 


325 


0455 


0750 


1595 


11.25/87 


329 


0438 


0734 


1854 



Table 6. OBSERVATION DAYS NOT USED IN THE ANALYSIS. 



Date 


Reason for not analyzing 


10/15/87 


No SEA-TAC weather observations 


10/16,87 


No satellites at Monterey for first 10 minutes 


10/23/87 


Slope distance a s > 10 m for first 10 minutes 


10/30/87 


Disk error 


10/31/87 


No SEA-TAC weather observations 


1 1/06/87 


Unhealthy satellites 


11/07/87 


Unhealthy satellites 


11/17/87 


Unhealthy satellites 


1 1/18/87 


Unhealthy satellites 


11/19/87 


Unhealthy satellites 


11/20/87 


Unhealthy satellites 


1 1/24/87 


No satellites at Monterey for first 20 minutes 
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Day 283 was processed with Monterey as the fixed reference station because the 
data set would not partition into 10-minute segments when Sand Point was used as the 
reference station. 

The days that were missing one or two segments were excluded from the analysis, 
so that changes in the sample variances would not be due to unequal sample populations 
between the segments. 

Reprocessing the first segment for 23 October failed to reduce the slope distance 
sigma to less than 10 m because that segment had only nine triple difference observa- 
tions because the receiver frequently lost lock on the satellites. That day was not used 
in the analysis, even though its other segments had slope distance sigma's less than 10 
m without reprocessing. 
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IV. RESULTS AND DISCUSSION 



A. GENERAL 

To study optimized times, five cases are studied: 

1 Process all data when four satellites at least 15° above horizon 

2 Begin processing when five satellites at least 15° above the horizon 

3 Begin processing as in 2, but delete fifth satellite 

4 Begin processing data 40 minutes later than in 1 

5 Begin processing data 70 minutes later than in 1 

Case 1 is essentially the processing of the full data set. Each of the other cases is a 
subset of Case 1. Trim640 allows the user to designate at what time within the full ob- 
serving period that the data loading should begin. For Case 3, Trim640's ability to flag 
data was used to exclude the fifth satellite (SV 12) from the processing. Each day of 
Table 5 was processed for each of the five cases. 

During the course of the discussions it will be necessary’ to distinguish between ob- 
servation periods and the time of the observations fixed with respect to satellite geom- 
etry. As the satellite geometry (PDOP) begins with the availability of four satellites, the 
time of observations can be defined in terms of the Case 1 start time. Observation pe- 
riods are determined from the start time for each case. Times of observations will given 
as equivalent Case 1 times and is obtained by adding the case observation length to the 
case's time offset from Case 1 (Case 2 and 3's offsets are 20 minutes; Case 4, 40 minutes 
and Case 5, 70 minutes). 

B. ACCURACY 

1. Slope Distance 

The slope distance errors for Cases 1, 2, 4, and 5 are presented in Table 7. Then 
results show that accuracy to better than 1.0 ppm is achieved for any observation period, 
but that there are differences among the cases and with changing observation periods 
within each case. 

Cases 1, 2, and 4 exhibit similar behavior as the observation period increases - 
they become less positive (or more negative) as the observation period increases until 
they reach a minimum, then they reverse their trends and become less negative. Positive 
error indicates that the measured baseline is longer than the standard values, so that the 
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measured baselines are exhibiting an accordion effect as they shorten then lengthen as 
more observations are included in the solutions. While the minimum occurs after dif- 
ferent observing periods (140, 110, 90 minutes for Cases 1, 2, and 4 respectively), they 
occur at about the same absolute time with respect to the Case 1 start time (140, 130, 
and 130 minutes for Cases 1,2, and 4 respectively). As the start limes occur later with 
each case, the errors for shorter observing periods become less positive and the ranges 
of the errors for each case decrease (the range of Case 1 errors is + 36 to -18.4 cm; + 7.S 
to -18.1 cm, for Case 2; and -10.3 to -15.1 cm, for Case 4). The error after the entire 
observation session decreases with later observing start times ( -14.8 cm. -9.1 cm, and 
-0.4 cm for Cases 1, 2, and 4 respectively). 

Case 5 behaves similarly to the previous cases except that its minimum occurs 
after only 10 minutes of observations and adding more observations causes the error to 
become positive. The error is largest after the entire observing session (38.4 cm). Cases 
1, 2, and 4 start their observations prior to the first PDOP peak that occurs at 60 min- 
utes (Figure 4) while Case 5 starts after the PDOP peak, so that starting the observa- 
tions close to the larger PDOP peak and including the PDOP peak observations can 
reduce the error and the required observation period. The effect of the PDOP peak upon 
the mean slope distance error is readily apparent in comparing Cases 4 and 5. Case 4 
remains negative without the early positive error, and Case 5 remains positive as it lacks 
most of the observations from about the first PDOP peak. 

The first PDOP peak at 60 minutes differs from the second PDOP peak in that 
it has a higher value, is symmetric, and occurs farther from the four satellite observation 
start time than the second peak is from the four satellite observation stop time, i.e., the 
second peak occurs with lower elevation angles than the first which implies larger 
tropospheric and ionospheric errors. 

The results for Case 3 (excluding the fifth satellite) slope distance errors are 
presented in Table 8. Case 3 was studied because the Case 2 results showed a decrease 
in the initial slope distance errors, and Case 3 was to study the effects of observing five 
satellites. The results show that using only four satellites makes a difference of only a 
few centimeters from the results using all available satellites. It should be noted that 
Case 3 uses only three satellites once SV 6 sets after 100 minutes of observations. The 
Case 3 minimum, -11.5 cm, was less negative than the Case 2 minimum of -18.1 cm. 



25 



Table 7. MEAN SLOPE DISTANCE ERROR FOR CASES 1, 2. 4. AND 5 



Observation 

Period 

[min| 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


36.0 


0.29 


7.8 


0.06 


-10.3 


-0.08 


-2.3 


-0.02 


20 


30.4 


0.25 


2.5 


0.20 


-10.9 


-0.09 


5.7 


0.05 


30 


21.7 


0.18 


-0.4 


-0.00 


-12.2 


-0.10 


19.2 


0.16 


40 


15.1 


0.12 


-3.3 


-0.03 


-13.0 


-0.1 1 


19.7 


0.16 


50 


9.8 


0.08 


-6.0 


-0.05 


-12.9 


-0.11 


17.7 


0.14 


60 


5.3 


0.04 


-8.6 


-0.07 


-10.9 


-0.09 


17.9 


0.15 


70 


0.8 


0.0 1 


-10.4 


-o.os 


-13.2 


-0.11 


20.6 


0.17 


SO 


-3.0 


-0.02 


-10.7 


-0.09 


-15.0 


-0.12 


23.2 


0.19 


90 


-6.0 


-0.05 


-14.2 


-0.12 


-15.1 


-0.12 


27.2 


0.22 


100 


-7.8 


-0.06 


-16.8 


-0.14 


-13.4 


-0.11 


31.7 


0.26 


110 


-12.1 


-0.10 


-18.1 


-0. 14 


-11.2 


-0.09 


38.4 


0.31 


120 


-15.9 


-0.13 


-17.2 


-0.14 


-8.0 


-0.07 


- 


- 


130 


-18.1 


-0.15 


-16.0 


-0.13 


-4.8 


-0.04 


- 


- 


140 


-18.4 


-0.15 


-14.0 


-0.11 


-0.4 


-0.00 


- 


- 


150 


-18.0 


-0.15 


-11.7 


-0.10 


- 


- 


- 


- 


160 


-17.0 


-0.14 


-9.1 


-0.07 


- 


- 


- 


- 


170 


-15.6 


-0.13 


- 


- 


- 


- 


- 


- 


175 


-14.8 


-0.12 


- 


- 


- 


- 


- 


- 



Table 8. CASE 3 MEAN ERRORS 



Observation 

Period 

[min] 


Slope 

Distance 


AX 


AY 


AZ 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


8.8 


0.07 


-287.0 


-6.97 


177.6 


2.48 


-21.7 


-0.24 


20 


3.7 


0.03 


-247.0 


-6.00 


155.5 


2.17 


-15.6 


-0.17 


30 


-0.4 


-0.00 


-214.3 


-5.21 


126.3 


1.76 


-1.9 


-0.02 


40 


-4.6 


-0.04 


-199.4 


-4.84 


110.4 


1.54 


9.5 


0.10 


50 


-8.5 


-0.07 


-170.6 


-4.15 


87.6 


1.22 


19.7 


0.22 


60 


-10.9 


-0.09 


-142.2 


-3.44 


64.9 


0.91 


27.9 


0.30 


70 


-11.5 


-0.10 


-115.7 


-2.82 


44.8 


0.63 


32.6 


0.36 


so 


-10.1 


-0.08 


-90.4 


-2.19 


24.5 


0.34 


35.1 


0.3S 


90 


-10.3 


-0.08 


-58.7 


-1.43 


7.2 


0.10 


34.7 


0.38 


100 


-10.0 


-0.08 


-38.5 


-0.95 


9.2 


0.13 


23.7 


0.26 


110 


-10.3 


-0.08 


-15.8 


-0.39 


3.9 


0.05 


18.0 


0.20 


120 


-10.2 


-0.08 


5.9 


0.15 


-1.0 


-0.01 


11.8 


0.13 


130 


-9.7 


-0.08 


28.7 


0.70 


- 7.3 


-0.10 


5.9 


0.06 


140 


-8.8 


-0.07 


52.9 


1.29 


-18.6 


-0.26 


2.5 


0.03 


150 


-7.5 


-0.06 


73.7 


1.80 


-29.8 


-0.42 


0.3 


0.00 


160 


-5.2 


-0.04 


97.2 


2.36 


-45.0 


-0.63 


-1.5 


-0.02 
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2. Baseline Components 

The AX and AY accuracies are better than 10.0 ppm for all observing periods 
while AZ accuracy is better than 1.0 ppm for all observing periods. The accuracy of the 
baseline components is expected to be less than the accuracy of the baseline because the 
baseline errors are mostly perpendicular to the baseline itself [Remondi, 1984], The 
baseline component results (Tables 9, 10, and 11) show that the AX and AY errors are 
greater than the baseline components while the AZ errors are about the same order of 
magnitude as the slope distance errors. The AX and AY errors are negatively correlated 
which is the result of the correlations of the triple differences and the semi-circular tracks 
of the satellites (Figure 5). Case 3 (Table 8) shows little difference from Case 2 in the 
AX error, and a more negative AY error is offset by a less negative AZ error. 

The smallest mean errors for the baseline components are found at various ob- 
serving periods. Zero mean error for all the baseline components is achieved with fewer 
observations as the observing start times occur closer to and before the larger PDOP 
peak. 



Table 9. MEAN AX ERROR FOR CASES 1, 2. 4, AND 5 



Observation 

Period 

(min) 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


-252.4 


-6.13 


-296.3 


-7.19 


-136.3 


-3.30 


29.6 


0.73 


20 


-255.4 


-6.21 


-256.0 


-6.22 


-136.4 


-3.30 


83.9 


2.04 


30 


-255.5 


-6.21 


-212.1 


-5.15 


-107.6 


-2.62 


132.9 


3.23 


40 


-238.0 


-5.78 


-194.8 


-4.74 


-80.1 


-1.94 


183.9 


4.47 


50 


-211.2 


-5.12 


-167.8 


-4.08 


-52.7 


-1.29 


191.7 


4.66 


60 


-197.0 


-4.79 


-137.4 


-3.33 


-24.8 


-0.61 


192.9 


4.69 


70 


-171.8 


-4.17 


-107.3 


-2.60 


21.9 


0.53 


193.4 


4.70 


80 


-145.1 


-3.52 


-76.8 


-1.87 


48.4 


1.17 


195.4 


4.75 


90 


-116.9 


-2.84 


-31.8 


-0.78 


67.7 


1.65 


201.7 


4.91 


100 


-87.3 


-2.11 


3.1 


0.08 


82.3 


1.99 


208.7 


5.08 


110 


-45.2 


-1.09 


28.4 


0.70 


96.2 


2.33 


221.0 


5.37 


120 


-10.7 


-0.27 


48.6 


1.18 


112.0 


2.72 


- 


- 


130 


15.2 


0.36 


65.7 


1.60 


124.8 


3.04 


- 


- 


140 


35.6 


0.87 


83.3 


2.02 


139.9 


3.40 


- 


- 


150 


53.6 


1.30 


97.7 


2.37 


- 


- 


- 


- 


160 


70.7 


1.72 


110.1 


2.67 


- 


- 


- 


- 


170 


84.3 


2.04 


- 


- 


- 


- 


- 


- 


175 


91.7 


2.24 


- 


- 


- 


- 


- 


- 
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Table 10. MEAN AY ERROR FOR CASES 1, 2, 4, AND 5 



Observation 

Period 

[min] 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


188.0 


2.63 


199.2 


2.78 


110.0 


1.54 


9.1 


0.13 


20 


181.2 


2.53 


181.4 


2.53 


106.3 


1.49 


-28.0 


-0.39 


30 


179.3 


2.50 


152.4 


2.13 


93.6 


1.31 


-68.6 


-0.96 


40 


171.2 


2.39 


139.3 


1.95 


79.8 


1.11 


-95.2 


-1.33 


50 


153.0 


2.14 


125.3 


1.75 


66.6 


0.93 


-92.6 


-1.29 


60 


141.0 


1.97 


109.5 


1.53 


52.2 


0.73 


-91.9 


-1.28 


70 


127.7 


1.78 


94.3 


1.32 


29.3 


0.41 


-89.7 


-1.25 


80 


114.0 


1.59 


79.3 


1.11 


21.6 


0.30 


-89.3 


-1.25 


90 


100.5 


1.40 


58.4 


0.82 


11.1 


0.16 


-95.6 


-1.33 


100 


86.4 


1.21 


46.4 


0.65 


4.0 


0.06 


-105.4 


-1.47 


110 


66.7 


0.93 


34.0 


0.47 


-3.3 


-0.05 


-121.0 


-1.69 


120 


55.8 


0.78 


23.5 


0.33 


-15. S 


-0.22 


- 


- 


130 


43.3 


0.60 


14.8 


0.21 


-27.7 


-0.39 


- 


- 


140 


33.9 


0.47 


2.9 


0.04 


-42.5 


-0.59 


- 


- 


150 


25.5 


0.36 


-9.0 


-0.13 


- 


- 


- 


- 


160 


14.4 


0.20 


-17.0 


-0.24 


- 


- 


- 


- 


170 


4.0 


0.06 


- 


- 


- 


- 


- 


- 


175 


-1.8 


-0.03 


- 


- 


- 


- 


- 


- 



Table 1 1. MEAN AZ ERROR FOR CASES 1, 2, 4, AND 5 



Observation 

Period 

[min] 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


-82.2 


-0.90 


-33.2 


-03.6 


-11.0 


-0.12 


-17.4 


-0.19 


20 


-68.0 


-0.75 


-30.2 


-0.33 


- 7.3 


-0.08 


-23.6 


-0.26 


30 


-54.8 


-0.60 


-23.4 


-0.36 


-8.5 


-0.09 


-32.0 


-0.36 


40 


-47.3 


-0.52 


-17.0 


-0.19 


-8.9 


-0.10 


-34.8 


-0.38 


50 


-38.1 


-0.42 


-14.6 


-0.16 


-11. 1 


-0.12 


-37.7 


-0.41 


60 


-28.9 


-0.32 


-12.3 


-0.13 


-15.1 


-0.17 


-39.1 


-0.42 


70 


-23.8 


-0.26 


-11.6 


-0.13 


-15.1 


-0.17 


-44.7 


-0.49 


so 


-20.0 


-0.22 


-13.1 


-0.14 


-18.5 


-0.20 


-49.5 


-0.54 


90 


-18.1 


-0.20 


-12.4 


-0.14 


-18.9 


-0.20 


-52.7 


-0.58 


100 


-17.9 


-0.20 


-15.1 


-0.17 


-22.2 


-0.24 


-54.2 


-0.59 


110 


-15.6 


-0.17 


-15.1 


-0.17 


-25.7 


-0.28 


-56.5 


-0.62 


120 


-17.5 


-0.19 


-17.2 


-0.19 


-27.3 


-0.30 




- 


130 


-16.4 


-0.18 


-19.8 


-0.22 


-28.0 


-0.31 




- 


140 


-17.9 


-0.20 


-21.0 


-0.23 


-29.3 


-0.32 




- 


150 


-19.9 


-0.22 


-21.3 


-0.23 


- 


- 




- 


160 


-20.3 


-0.22 


-21.9 


-0.24 


- 


- 




- 


170 


-20.1 


-0.22 


- 


- 


- 


- 




- 


175 


-20.1 


-0.22 


- 


- 


- 


- 




- 
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NORTH 




Figure 5. Sky Plots of satellite tracks for Monterey: Elevation angles are dotted 

concentric circles. Zenith is at the center. 
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C. REPEATABILITY 
1. Slope Distance 

The day-to-day repeatabilities, represented as the 2a level, for Cases 1, 2, 4, and 
5 are presented in Table 12, and for Case 3, in Table 13. All the cases achieve 1.0 ppm 
repeatability for any observing period except Case 1 which requires 20 minutes of ob- 
servations. Repeatability eventually reaches better than 0.5 ppm after 60 minutes of 
observations for any case. 

The minimum 2cr levels for all cases are reached at the 80 to 90 minute time of 
observation, which is about 30 minutes after the larger PDOP peak. A slight increase 
in the 2a level is centered about the 120 to 130 minute time of observation for all cases 
which is near the PDOP minimum, after which the 2a level decreases slightly as obser- 
vations from the second PDOP peak are included in the solutions. Case 4 had the 
narrowest range of 2 a values, and Case 1 had the widest range of 2 a values. 



Table 12. SLOPE DISTANCE 2-SIGMA VALUES FOR CASES 1, 2, 4, AND 5 



Observation 

Period 

[min] 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


176.8 


1.44 


64.8 


0.53 


57.8 


0.47 


52.0 


0.42 


20 


118.2 


0.96 


63.2 


0.5! 


56.8 


0.46 


53.4 


0.43 


30 


84.1 


0.68 


58.6 


0.4S 


59.2 


0.48 


60.8 


0.49 


40 


69.8 


0.57 


58.4 


0.47 


53.8 


0.44 


57.0 


0.46 


50 


64.0 


0.52 


58.6 


0.48 


54.4 


0.44 


55.4 


0.45 


60 


61.8 


0.50 


55.2 


0.45 


57.6 


0.47 


55.0 


0.45 


70 


60.4 


0.49 


55.0 


0.44 


59.6 


0.48 


53.6 


0.44 


80 


57.4 


0.47 


56.8 


0.46 


60.4 


0.49 


53.6 


0.44 


90 


56.2 


0.46 


59.0 


0.48 


60.8 


0.49 


57.6 


0.46 


100 


57.2 


0.47 


60.0 


0.49 


59.6 


0.48 


62.6 


0.51 


110 


59.2 


0.48 


60.4 


0.49 


57.8 


0.47 


67.8 


0.55 


120 


60.0 


0.49 


59.6 


0.48 


56.2 


0.46 


- 


- 


1 30 


61.0 


0.50 


59.0 


0.48 


56.0 


0.46 


- 


- 


140 


60.6 


0.49 


58.2 


0.47 


55.4 


0.45 


- 


- 


150 


60.0 


0.49 


57.8 


- 


. 


- 


- 


- 


160 


59.0 


0.48 


58.4 


- 


- 


- 


- 


- 


170 


58.6 


0.48 


- 




. 


- 


- 


. 


175 


58.6 


0.48 


- 




- 


- 


- 


- 
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Table 13. CASE 3: 2-SIGMA VALUES 



Observation 

Period 

[min] 


Slope 

Distance 


AX 


AY 


AZ 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


62.4 


0.51 


326.0 


7.92 


140.8 


1.97 


76.6 


0.84 


20 


59.4 


0.49 


304.2 


7.39 


135.8 


1.90 


80.8 


0.89 


30 


57.8 


0.47 


3 IS. 8 


7.75 


109.0 


2.36 


76.6 


0.S4 


40 


58.0 


0.47 


301.0 


7.31 


173.4 


2.42 


S2.2 


0.90 


50 


58.8 


0.48 


287.4 


6.98 


195.S 


2.60 


95.4 


1.04 


60 


54.0 


0.44 


272.2 


6.61 


197.0 


2.75 


97.0 


1.06 


70 


53.0 


0.43 


251.8 


6.12 


199.8 


2.79 


95.8 


1.05 


SO 


53.6 


0.44 


244.6 


5.94 


199.4 


2.79 


89.0 


0.98 


90 


56.4 


0.46 


235.6 


5.72 


187.0 


2.61 


83.4 


0.91 


100 


58.0 


0.47 


219.8 


5.34 


165.8 


2.32 


76.2 


0.S4 


110 


59.2 


0.48 


199.2 


4.84 


138.8 


1.94 


71.8 


0.79 


120 


59.0 


0.4S 


18S.8 


4.59 


127.4 


1.78 


6S.2 


0.75 


130 


59.0 


0.4S 


184.2 


4.48 


122.6 


1.71 


65.0 


0.71 


140 


58.8 


0.48 


1S8.4 


4.58 


127.2 


1.7S 


65.4 


0.72 


150 


58.8 


0.48 


190.4 


4.62 


128.4 


1.79 


64.8 


0.71 


160 


57.8 


0.47 


203.0 


4.93 


137.0 


1.91 


63.8 


0.70 



2. Baseline Components 

All the baseline components have repeatabilities better than 10.0 ppm for any 
observing period except for the Case 1 AX, which required 20 minutes of observations 
(Tables 13, 14, 15, and 16). 

It is interesting to note that while the AX 2 a values for the first segment of 
Cases 2, 3, 4, and 5 are less than Case I's first segment, the final Case 1 2 a value is less 
than the final segment of any other Case. The final 2 a values for Cases 2, 3, 4, and 5 
are greater than the Case 1 2 a values after an equivalent number of observations. 

The minimum la levels for AX and AY occur at or near the end of the observing 
sessions for Cases 1,2, 4 and 5. For Case 3 the minimum 2 a levels occur after 130 min- 
utes of observations. The minimum AZ 2 a level occurs after various observation peri- 
ods, but generally in the vicinity of 130 to 150 minute observation time, which is between 
the PDOP minimum and the second PDOP peak. Case 3 behaves in an opposite fashion 
from the other cases in that its minimum AX and AY 2 a levels occur at the 150-minute 
time of observation while its AZ minimum occurs at the end of the observation period. 



31 



Table 14. AX 2-SIGMA VALUES FOR CASES 1, 2, 4, AND 5 



Observation 

Period 

[min] 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


646.4 


15.69 


347.2 


8.45 


422.6 


10.27 


395.2 


9.60 


20 


531.4 


12.93 


311.4 


7.57 


343.4 


8.34 


375.4 


9.12 


30 


345.0 


8.41 


312.2 


7.58 


312.2 


7.32 


375.0 


8.87 


40 


248.2 


6.03 


311.0 


7.31 


286.8 


6.97 


336.2 


8.17 


50 


239.6 


5.83 


286.0 


6.95 


251.0 


6.10 


315.8 


7.43 


60 


230.4 


5.59 


269.2 


6.54 


246.8 


6.00 


272.8 


6.63 


70 


221.2 


5.39 


246.6 


5.99 


233.4 


5.67 


252.8 


6.14 


80 


209.2 


5.10 


243.2 


5.91 


212.8 


5.17 


231.0 


5.61 


90 


195.0 


4.74 


221.4 


5.38 


188.4 


4.43 


224.4 


5.45 


100 


202.0 


4.91 


194.8 


4.73 


16S.8 


4.10 


206.6 


5.02 


110 


180.8 


4.39 


168.8 


4.10 


159.2 


3.87 


203.4 


4.94 


120 


162.4 


3.94 


155.2 


3.77 


155.6 


3.78 


- 


- 


130 


141.4 


3.44 


146.4 


3.56 


141.4 


3.43 


- 


- 


140 


134.6 


3.26 


144.4 


3.51 


139.8 


3.40 


- 


- 


150 


129.0 


3.11 


133.0 


3.23 


- 


- 


- 


- 


160 


129.4 


3.14 


134.4 


3.27 


- 


- 


- 


- 


170 


121.8 


2.96 


- 


- 


- 


- 


- 


- 


175 


121.8 


2.96 


- 


- 


- 


- 


- 


- 



Table 15. AY 2-SIGMA VALUES FOR CASES 1, 2, 4, AND 5 



Observation 

Period 

[min] 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


110.8 


1.55 


204.0 


2.85 


244.2 


3.41 


325.0 


4.54 


20 


130.4 


1.82 


191.2 


2.67 


224.2 


3.13 


311.4 


4.34 


30 


140.8 


1.97 


190.6 


2.67 


220.4 


3.08 


298.4 


4.17 


40 


130.8 


1.83 


193.2 


2.70 


221.0 


3.09 


264.8 


3.70 


50 


135.4 


1.89 


192.2 


2.69 


212.0 


2.96 


237.0 


3.31 


60 


142.8 


1.99 


188.6 


2.63 


203.8 


3.00 


200.0 


2.79 


70 


144.4 


2.02 


185.2 


2.59 


214.6 


2.85 


181.4 


2.53 


80 


144.6 


2.02 


185.0 


2.58 


181.6 


2.54 


170.2 


2.38 


90 


145.4 


2.03 


173.4 


2.42 


147.4 


2.06 


177.2 


2.48 


100 


145.0 


2.10 


152.4 


2.13 


133.8 


1.87 


174.6 


2.44 


110 


139.0 


1.94 


127.6 


1.78 


124.6 


1.74 


182.2 


2.55 


120 


126.8 


1.77 


115.2 


1.61 


124.8 


1.74 


- 




130 


108.0 


1.51 


110.6 


1.55 


119.0 


1.66 


- 




140 


102.6 


1.43 


112.6 


1.57 


120.4 


1.68 


- 




150 


100.0 


1.40 


109.0 


1.52 


- 


- 


- 




160 


101.8 


1.42 


118.0 


1.65 


- 


- 


- 




170 


100.0 


1.40 


- 


- 


- 


- 


- 




175 


101.8 


1.42 


- 




- 


- 


- 





32 



Table 16. AZ 2-SIGMA VALUES FOR CASES 1, 2, 4, AND 5 



Observation 

Period 

[min] 


Case 1 


Case 2 


Case 4 


Case 5 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


73.6 


0.81 


65.2 


0.76 


92.8 


1.02 


88.0 


0.97 


20 


7S.0 


0.S6 


73.4 


0.81 


87.4 


0.96 


74.6 


0.S2 


30 


53.8 


0.59 


69.8 


0.77 


82.8 


0.91 


65.2 


0.72 


40 


54.6 


0.60 


71.6 


0.79 


78.4 


0.86 


57.8 


0.63 


50 


55.4 


0.61 


71.6 


0.79 


72.0 


0.79 


49.2 


0.53 


60 


60.2 


0.66 


68.4 


0.75 


68.4 


0.75 


47.8 


0.52 


70 


62.6 


0.69 


63.8 


0.70 


65.0 


0.71 


52.2 


0.57 


SO 


61.0 


0.67 


59. S 


0.70 


59.8 


0.66 


53.2 


0.58 


90 


59.2 


0.65 


58.2 


0.64 


55.2 


0.61 


56.2 


0.62 


100 


58.0 


0.64 


55.4 


0.61 


53.2 


0.58 


58.2 


0.64 


110 


56.6 


0.62 


53.4 


0.59 


52.2 


0.57 


59.4 


0.65 


120 


55.8 


0.61 


52.8 


0.58 


53.6 


0.59 


. 


- 


130 


54.4 


0.60 


52.4 


0.58 


54.0 


0.59 


- 


. 


140 


53.8 


0.59 


54.2 


0.59 


53.6 


0.59 


- 


_ 


150 


53.6 


0.59 


54.4 


0.60 


_ 


- 


. 


. 


160 


56.0 


0.61 


53.4 


0.59 


- 


- 


- 




170 


56.6 


0.62 


- 


- 


- 


- 


• 


- 


175 


56.0 


0.61 


- 


- 


- 


- 


* 


- 



3. Standard Deviation of the Mean 

The repeatability values can be used to estimate the standard deviations of the 
mean errors given in Tables 7 through 11 by using Equation (5.2). The values of Tables 
12 through 16 should be divided by v 28 (where 28 is the sample population) to compute 
the standard deviations of the means (at the la level). Generally, the repeatabilities were 
about five times the magnitudes of the mean errors; therefore, the uncertainties of the 
mean errors are on the order of the mean errors themselves. Allowing for the 0. 1 ppm 
uncertainty in the baseline and the baseline components and for the possible standard 
deviation of the means, accuracies to better than 1.0 ppm for the slope distances and 
10.0 ppm for the baseline components remain valid. 

D. ERROR EFFECTS 
1. 7- Day Means 

Because of the observations were made over a long period of time, Case 1 was 
subdivided into four groups comprised of seven consecutive observation days to study 
trends in the slope distance error to identify the contribution of various error sources to 
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optimizing observing times. The results for the slope distance error are presented in 
Table 17. 



Table 17. CASE 1 ERROR: 7-DAY MEANS 



Observation 

Period 

[min| 


Group 1 


Group 2 


Group 3 


Group 4 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


SO. 5 


0.65 


-19.0 


-0.15 


53.5 


0.43 


29.0 


0.24 


20 


62.6 


0.51 


21.1 


0.17 


26.9 


0.22 


10.9 


0.09 


30 


42.8 


0.35 


1.4 


0.01 


26.3 


0.21 


16.4 


0.13 


40 


30.6 


0.25 


-8.0 


-0.07 


27.3 


0.22 


10.4 


0.08 


50 


24.3 


0.20 


-12.9 


-0.10 


22.0 


0.18 


6.1 


0.05 


60 


21.5 


0.17 


-18.2 


-0.15 


16.1 


0.13 


1.8 


0.01 


70 


17.8 


0.14 


-24.6 


-0.20 


10.7 


0.09 


-0.8 


-0.01 


80 


13.2 


0.11 


-25.4 


-0.21 


6.1 


0.05 


-5.9 


-0.05 


90 


10.2 


0.08 


-27.6 


-0.22 


2.6 


0.02 


-9.2 


-0.07 


100 


10.7 


0.09 


-30.9 


-0.25 


0.8 


0.01 


-12.0 


-0.10 


110 


6.4 


0.05 


-36.4 


-0.30 


-2.4 


-0.02 


-16.2 


-0.13 


120 


2.0 


0.02 


-40.7 


-0.33 


-6.1 


-0.05 


-19.0 


-0.15 


130 


-0.6 


-0.00 


-43.6 


-0.35 


-6.6 


-0.05 


-21.6 


-0.18 


140 


-2.1 


-0.02 


-43.5 


-0.35 


-6.3 


-0.05 


-21.8 


-0.18 


150 


-2.3 


-0.02 


-42.5 


-0.35 


-5.6 


-0.05 


-21.8 


-0.18 


160 


-1.6 


-0.01 


-40.6 


-0.33 


-4.8 


-0.04 


-21.0 


-0.17 


170 


-0.2 


-0.00 


-38.3 


-0.31 


-3.4 


-0.03 


-20.5 


-0.17 


175 


0.7 


0.01 


-37.4 


-0.30 


-2.3 


-0.02 


-20.1 


-0.16 



The seven-day mean slope distance errors remain below 1.0 ppm for all groups 
for any observation period. The differences between the groups are mainly in the 
predominance of negative errors in Groups 2 and 4 while Groups 1 and 4 have pre- 
dominantly positive errors. 

To account for the change in the characters of the groups, the change in the 
major sources of error are examined. 

2. Ephemeris Errors 

Examining the AODE (Age of Data Ephemeris) for all satellites during the 
eight-week observing period, showed that the oldest AODE for a single satellite was nine 
hours, but most were less than five hours. Averaging the AODE for all satellites for 
each observing day showed a range of seven hours during the eight weeks of observa- 
tions (Figure 6). Assuming 12 hours is the largest AODE, the maximum ephemeris 
linearization error would be 50 cm. Using Equation (2.3) results in 0.025 ppm baseline 
relative error. 
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Figure 6. Mean age of data (ephemeris) 

The changes in the seven-day mean AODE's do not correspond well to the 
changes in the groupings, especially as the largest change in AODE from group 3 to 4 
does not correspond to a similar change in the group 3 to 4 mean slope distance errors. 
This is not surprising because while the change is relatively large, the magnitude of the 
orbital errors is expected to be small. 

The ephemeris error will appear as a bias during any one observing session. By 
averaging over eight weeks some of the errors will cancel and appear as variance. Be- 
cause of the short observing sessions, little variation in the ephemeris error is expected 
during any single sessions 

3. Ionospheric Errors 

All the observing sessions completed well before 1600 local standard time and 
approached the normally ionospheric activity minimum of 0600 (Figure 7). The near 
north-south orientation of the Monterey-Sand Point baseline places both stations in the 
same time zone, so that ionospheric errors between the two stations will be correlated 
to some degree and reduced in the differencing. 
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Figure 7. Difference between observation end time and 0600 PST 




Figures. Electron fluenee: Units electrons — cm~ 2 —day' —sr~'. 
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The weekly Preliminary Report and Forecast of Solar Geophysical Data | Space 
Environment Services Center, 1987] for the eight weeks of the observations described the 
solar and geomagnetic activity as generally quiet or low. 

Though the observing period occurred during a general lull in ionospheric ac- 
tivity, there are two ionospheric cfl'ects that can be easily examined: the diurnal effect 
as the observing periods occur closer to the diurnal ionospheric activity minimum and 
sporadic effects from electron fluxes. The change in ionospheric range error for the 
eight-week observing period was computed from Henson and Collier [1986, Equation 
(3)j using the difference in Total Electron Content (TEC) determined from Spilker [1980, 
Figure 1-11] and Henson and Collier [1986, Figure 2] for the mean observation times for 
the first and last groups. The maximum change in range error, because of diurnal 
ionospheric changes from the first to the last week, is approximately 1.5 m which by 
Equation (2.3) is equivalent to a change in the relative baseline error of only 0.075 ppm. 

While the 0.075 ppm trend is not discernible because of the 0. 1 ppm uncertainty 
in the baseline and the larger uncertainties in the measurements, the mean number of 
observations (Table 5) increased with each later group. The mean number of triple 
difference observations available were : 1519 ±45, 1522 ±22, 1601 ±124, and 1726 ±174 
for the first through the last groups respectively. The increased number of observations 
could be attributed to stronger signal-to-noise ratio with the reduced ionospheric 
dispersion as the observations approach the diurnal ionospheric activity minimum. 

The group to group changes did not show an apparent trend, but alternated 
between predominantly negative and positive mean slope distance errors indicating a 
sporadic effect, or complex interaction between the error sources. To examine sporadic 
ionospheric effects, daily Electron Fluence (EF) from the weekly Preliminary Report and 
Forecast of Solar Geophysical Data (Space Environmental Services Center, 1987], and 
seven-day means of the EF were computed (Figure 8). Electron fluence 
( electrons — cnr 2 —day- 1 — sr _1 ) is the daily average of electron flux with energies 
greater than 2 Mev as measured by the GOES-7 satellite. High values of the seven-day 
mean EF corresponded to predominantly positive seven-day mean slope distance errors 
(Groups 1 and 3), and low mean EF corresponded to negative mean slope distance error. 

4. Tropospheric Errors 

Refractivity, A r , was computed using the meteorological parameters from the 
time closest to the middle of each observation session and Remondi [1984, equations 
2.30 and 2.31]. Possible indicators of the influence of the troposphere upon the baseline 
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solution are the difference in refractivity, AA\ between the baseline stations (Figure 9) 
and the mean A’ of both stations (Figure 10). The difference can indicate the level of 
correlation of the tropospheric errors. Correlated tropospheric errors are reduced in 
differencing. The mean A’ (A') will indicate the magnitude of the refraction of the carrier 
signal to both sites. 

The changes in the difference in the seven-day mean A’ (Figure 9) did not have 
an apparent correspondence to the changes in the seven-day group mean slope distance 
errors. Instead A' shows a pattern similar to the group slope distance errors 
(Figure 10). High mean refractivity corresponded to positive mean slope distance errors 
(Groups 1 and 3), and low mean refractivity corresponded with negative mean slope 
distance errors. The difference between the high and the low mean A' are small and do 
not readily account for the magnitude of the changes. 




Figure 9. Difference in refractivity between Monterey and Sand Point 
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Figure 10. Mean Monterey-Sand Point refractivity 
5. 7-Day Repeatability 

The Group repeatabilities could not be given a general characterization as the 
Group means were because the was no consistent pattern to the Group-to-Group 
changes in their repeatabilities. For completeness, the la values for the various error 
sources are presented in Table 18, and the Group la values are presented in Table 19. 



Table 18. GROUP 2-SIGMA VALUES FOR THE ERROR SOURCES 



Error Source 


Group 1 


Group 2 


Group 3 


Group 4 


AODE 

[hours] 


1.6 


2.0 


1.6 


2.4 


Observation End Time 
- 0600 PST 
[minutes] 


24.8 


45.0 


37.6 


56.8 


Electron Fluence 
[electrons — enr 2 —day' 1 — sr -1 ] 


610.0 


59.6 


294.8 


21.4 


AN 


27.4 


21.2 


27.8 


23.8 


.V 


15.2 


10.6 


9.0 


17.8 
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Table 19. CASE 1: 2-SIGMA VALUES, 7-DAY GROUPINGS 



Observation 

Period 

[niin] 


Group 1 


Group 2 


Group 3 


Group 4 


cm 


ppm 


cm 


ppm 


cm 


ppm 


cm 


ppm 


10 


129.8 


1.06 


81.2 


0.66 


155.4 


1.26 


232.0 


1.89 


20 


S9.4 


0.73 


110.2 


0.90 


141.8 


1.15 


124.8 


1.01 


30 


61.8 


0.50 


64.0 


0.52 


108.2 


0.88 


89.6 


0.73 


40 


60.8 


0.49 


56.2 


0.46 


80.6 


0.70 


64. S 


0.53 


50 


58.0 


0.47 


53.6 


0.44 


72.6 


0.59 


52.4 


0.43 


60 


59.8 


0.49 


53.4 


0.43 


65.4 


0.53 


45.6 


0.37 


70 


62.0 


0.50 


46.2 


0.38 


60.2 


0.49 


44.8 


0.36 


SO 


57.0 


0.46 


49.0 


0.40 


53.4 


0.43 


48.2 


0.39 


90 


57.0 


0.46 


47.4 


0.39 


51.4 


0.42 


47.6 


0.39 


100 


55.6 


0.45 


46.0 


0.37 


50.4 


0.41 


50.0 


0.41 


110 


57.4 


0.47 


46.8 


0.38 


50.4 


0.41 


53.2 


0.43 


120 


61.8 


0.50 


46.8 


0.38 


48.2 


0.39 


55.4 


0.45 


130 


64.6 


0.53 


44.4 


0.36 


45.4 


0.37 


58.0 


0.47 


140 


65.0 


0.53 


43.0 


0.35 


43.4 


0.35 


60.2 


0.49 


150 


64.4 


0.53 


43.4 


0.35 


42.2 


0.34 


60.2 


0.49 


160 


63.4 


0.52 


42.8 


0.35 


42.6 


0.35 


59.2 


0.48 


170 


62.8 


0.51 


42.6 


0.35 


43.8 


0.36 


5S.6 


0.48 


175 


63.0 


0.51 


41.6 


0.34 


45.0 


0.37 


57.6 


0.47 



E. EFFECTS OF THE C/A CODE 

The best C/A code position for Monterey was extracted from each observing session 
for comparison with the carrier phase results. The coordinates of Sand Point were sub- 
tracted from the 28-day mean of the Monterey C/A code coordinates to form the 
baseline components. The slope distances were computed from the square root of the 
sum of the squares of the baseline components. The mean errors of the C/A code de- 
rived baseline were on the order of the carrier phase mean errors, but the 7a levels are 
several times greater than the carrier phase results (Table 20). These results are for one 
instantaneous position determination and would improve with time averaging within the 
observing period or Doppler-smoothing. 

The C/A code has its most significant impact upon the differenced carrier phase 
solutions when there are few carrier phase observations. Table 21 shows that reproc- 
essing the carrier phase observations was required for several segments because either 
more than ten percent of the observations were rejected or the slope distance solution 
sigma was greater than 10 m, and that all the reprocessed segments were confined to the 
first two segments. 
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Table 20, BEST C/A CODE RESULTS AND ERRORS 



Component 


Distance (rn) 


Relative Error (ppm) 


AX 


-411587.334 


-8.2 -30.0 


AY 


-715776.997 


. 0.8 +2S.5 


AZ 


-91 1739.358 


-3.3 +17.7 


Slope Distance 


1230044.511 


-0.6 ±7.2 



Table 21. CASE 1: DATA SEGMENTS BEFORE AND AFTER REPROCESS- 
ING 



Day 


Segment 


Total 

Obs 


Before 


After 


Rejected 

Obs 


Slope 

Distance 

(m) 


Rejected 

Obs 


Slope 

Distance 

(m) 


Edit 

Multiplier 


Case 1 


272 


1 


62 


14 


1230050.192 


0 


1230045.846 


99 


275 


1 


42 


17 


1230268. 9S8 


0 


1230045.294 


5 


276 


1 


51 


5 


1230047.276 


0 


1230047.020 


99 


280 


1 


61 


41 


1230227.492 


0 


1230045.717 


99 


280 


2 


110 


44 


1230093.548 


0 


1230045.915 


5 


287 


1 


48 


6 


1230048.636 


0 


1230045.326 


99 


294 


■> 

Am 


124 


23 


1230044.230 


12 


1230046.146 


5 


308 


1 


45 


* 


1230042.920 


0 


1230045.972 


99 


Case 2 


314 


1 


78 


56 


1230074.209 


0 


1230045.282 


3.5 


Case 3 


314 


1 


59 


37 


1230074.209 


0 


1230045.373 


3.5 


Case 4 


281 


1 


78 


22 


1230034.198 


0 


1230045.082 


3.5 


281 


2 


169 


47 


1230041.504 


11 


1230045.031 


3.5 


282 


1 


62 


20 


1230039.213 


0 


1230045.024 


3.5 


Case 5 


287 


2 


159 


43 


1230043.418 


3 


1230045.177 


3.5 


308 


1 


86 


24 


1230043.047 


0 


1230045.472 


3.5 


314 


1 


52 


* 1 


1230051.620 


0 


1230048.000 


99.0 



The carrier phase solutions were improved by the least-squares processing of the 
pseudoranges to obtain better initial antenna position estimates and also to correct the 
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carrier phase time-tags. All the Case 1 reprocessed segments were improved by forcing 
the processing program to accept more carrier phase observations by increasing the edit 
multiplier. 

Because of the success of improving the solutions by the pseudorange pre- 
processing, all segments for three days (276, 287, and 294) were reprocessed to see if their 
already acceptable solutions would improve. The solutions either did not change or 
became slightly less accurate. 

F. COMPARISON WITH PREVIOUS STUDIES 

Remondi [1984] concluded that 100 ppm accuracy was achievable in 30 minutes and 
10 ppm, in one hour - regardless of baseline length - based upon single-frequency, single 
difference solutions using precise ephemerides over baselines less than 100 km. To study 
the required observation times for various accuracies, he partitioned his data into 
15-minute observation spans. While decorrelated triple differences would be capable of 
performing relative geodesy at 1 ppm level, correlated triple differences may provide 5-10 
ppm or better [Remondi, 1985]. 

Some long baseline surveys include: 

Bock et al. [19S4] achieved 2 ppm accuracy in closure of a transcontinental net aided 
by external atomic clocks using 10 hours of single frequency, single difference obser- 
vations. 

Cannon et al. [1985] achieved 0.1 to 0.7 ppm repeatability using single frequency 
non-differenced network solutions without ionospheric correction for 1700-km 
baselines between two California stations and Calgary, Alberta, Canada over a four- 
day period. The range in repeatability was attributed to satellite clock errors. They 
also found that the largest errors occurred in the AX component. 

Mader and Abell [1985] found that single frequency long baseline GPS results using 
single differenced carrier phase observables agreed to one ppm with VLBI measured 
baseline lengths. Their 2-day repeatability was 0.24 ppm. 

Goad et al [1985] measured a 1302-km baseline between California and Texas to better 
than 1.0 ppm compared to the VLBI measured distance, and their two-day repeat- 
ability was better than 0.5 ppm. 

Bertiger and Lichten [1987] achieved 6 parts per billion repeatability over a 1314-km 
east-west baseline over separate four and seven-day periods using dual-frequency re- 
ceivers with a fiducial network for orbit determination, and some water vapor 
radiometers. 

Except for the Bertiger and Lichten [1987] results, the Case 1 results are competitive 
with the above studies despite the inherent low resolution and the added burden of the 
correlations of the correlated triple differences. A possibility for the optimum results of 
Case 1 is that a good solution is insured by setting both receivers to their known WGS84 
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coordinates. For real-time surveys, one end of the baseline will usually not be known 
to one meter accuracy. 

During the eight-week observation period, the third 4000SX was delivered to NFS. 
The antenna was installed on the same rooftop, but approximately 35 m to the north 
of the Monterey antenna. Two days of data using Monterey as- the reference station 
were used to fix the location of the new antenna. On Day 294, the new receiver's coor- 
dinates were offset 37 m, equivalent to 30 ppm baseline error. The difference between 
the offset antenna and Monterey solutions for the baseline to Sand Point is 0.7 ppm for 
the first segment (Figure 11). 




Figure 11. Results of station 2 offset 

For the remaining segments the difference is about 0.05 ppm. The first segment will 
be more affected by receiver offset because the small number of carrier phase observa- 
tions will be sensitive to the initial position estimate provided by the C/A code. The C/A 
code, in turn, will be sensitive to the receiver coordinates because the receiver will com- 
pute the slant range to the satellites using those coordinates. 

G. MISSING EPOCHS 

Missing epochs prevented a sixth case study in which the number of observations 
are reduced by flagging the data. The sixth case study was to have studied the effects 



43 




Figure 12. Fraction of available triple difference observations 

of reducing the number of observations without changing the geometry. Several days 
of data were processed in which only every other epoch was to be used in the processing. 
An hour of data was needed to obtain enough triple difference observations to consist- 
ently have less than 10% of the triple difference observations rejected and slope distance 
o ' s less than 10 m. This inability to form a reliable solution indicated a lack of carrier 
phase observations, so the number of observations were counted. 

The mean number of observations (accumulative) available for each segment, as a 
fraction of the maximum number of triple difference observations remained steady until 
120 minutes when it increased (Figure 12). The maximum number of triple difference 
observations is: \l TD = («, — l)(« r —1) where: 

M td maximum number of triple difference observations 
n, number of observations epochs 

n s number of observed satellites 

Solutions using only 10 minutes of data at a time were examined to investigate the 
rise in available observations after 120 minutes. A curious phenomenon was discovered 
that while SV 6 was being tracked at Sand Point between five and ten of every 40 epochs 
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were lost (Figure 13). As soon as SV 6 was no longer being tracked at Sand Point, no 
complete loss of epochs was observed (Figure 14). 



START DATE/TIME: 1987/10/21 16:43:07. DAY OF YEAR 294 TOW 319360. 

STOP OATE/TIME: 1987/10/21 16:52:60. DAY OF YEAR 294 TON 319980. 

DATA AVAILABLE 
STATION: SAND POINT 

SV 6 I .. | 

SV 9 | | 

SV11 | | 

SV12 | | 

SV13 I | 

STATION: MONTEREY 

SV 6 | | 

SV 9 | | 

SV11 | | 

SV12 | | 

SV13 I | 



Figure 13. Missing epochs while tracking SV 6: Each dot represents one obser- 

vation: each column, one epoch. Sand Point is missing epochs 4. 8, 14, 
18, 26, 30, 34, and 40. From Trim640 output. 



START OATE/TIME: 1987/10/21 17:03:04. DAY OF YEAR 294 TOW 320580. 

STOP DATE/TIME: 1987/10/21 17:12:60. 0AY OF YEAR 294 TOW 321180. 

0ATA AVAILABLE 
STATION: SAN0 POINT 

SV 9 | I 

SV11 | I 

SV12 | I 

SV13 | I 

SV 6 | I 

STATION: MONTEREY 

SV 9 | I 

SVH | I 

SV12 I I 

SV13 I I 

SV 6 I I 



Figure 14. Continuous tracking without SV 6.: No missing epochs at Sand Point. 



The phenomenon is not multipath interference because it affects all satellite at each 
epoch. A simultaneous outage would suggest interference from some other emitter at 
or near 1540 MHz, if the outages continued throughout the observation session rather 
than being limited to when SV 6 was being tracked. No physical cause for the outages 
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is readily available, and the data have been forwarded to Trimble Navigation Ltd for 
further evaluation. 

H. DILUTION OF PRECISION AND RANGE ERRORS 

The case studies showed that the accuracy and precision of the slope distance were 
improved by using observations associated with the infinite PDOP peak rather than ei- 
ther the PDOP minimum or the secondary PDOP rise near the end of the observing 
sessions, and that improvements in the slope distance accuracy could be achieved using 
fewer observations when observations are made through the infinite PDOP peak. Pre- 
cision and accuracy of the solutions are dependent on geometry (fixed with respect to 
Case 1 time) and error sources rather than the number of observations (fixed with re- 
spect to length of observing period). 

The mean slope distance errors were governed by the ionospheric and tropospheric 
errors opposing each other. Early and late observations had positive mean errors while 
the errors about the infinite PDOP peak were negative. Tropospheric errors are partially 
corrected by the modified Hopfield model while the ionospheric errors remain uncor- 
rected, and as the satellites reach the infinite PDOP peak the full impact of the 
ionospheric errors project upon the baseline. 

The role reversal of the PDOP peak from point positioning to relative positioning 
requires some investigation. Because PDOP reflects the strength of the satellite geom- 
etry for a single station, the investigation should begin with the satellite-baseline (two 
station) geometry by examining a simple projection of errors upon a baseline using a 
single satellite whose track is directly over the baseline (Figure 15). 

From the law of cosines: 



b 



T~ Pi +p\~ lp\P7 cosy 



(6.1) 



where: 

b T true baseline 

p, range from station 1 to the satellite 

p 2 range from station 2 to the satellite 

y angle between p, and p 2 
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Similarly we define a measured baseline, b M : 

bit = (Pi + E P \f + (P 2 + E P i) 2 ~ 2(Pi + £J(p 2 + £ p2 ) cos y (6.2) 



where: 

error in pj 
E a error in p, 



Subtracting the square root of (6.1) from the square root of (6.2), and letting 
cosy = 1, gives: 



h± ~ b r / 1 , l2 j 

br V bl b\ 



(6.3) 



where A E, = E f 2 — E fl and A p = p ; — p, . From (6.3), Ap acts as an amplification factor 
for the difference in the range errors ( A E f ). As a satellite approaches the mid-point 
of a baseline, Ap -» 0 and b T >> A E„. 
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Similar arguments can be extended to point positioning using two satellites at one 
epoch or two epochs of one satellite by inverting the geometry of Figure 15, but the 
reduction in range errors would apply to point positioning as well. 

The effects of y have been ignored because y is small. From the law of sines : 



sin y = 



b sin 9 

P 



(6.4) 



where 0 = e + v. e is the elevation angle to the satellite and v is the vertical angle be- 
tween the slope distance and the local horizon. Expanding sin 9 as 
sin e cos v + cos e sin v , and ( sin e cos v) > > ( cos e sin v) because cos leal ( v is -5° for 
the 1230 km baseline ), then sin 9 can be approximated by sine. Approximating b by 
1230 km and p by 20000 km. then, by the law of sines, y is approximately 1° when e is 
15°; by the law of cosines, y is approximately 2.5° when the satellite is at the baseline 
midpoint. 

With cos 1°~0.9998 and cos 2.5 o ^0.9991, a very insignificant change, most of the 
deviation of the measured baseline from the true baseline must come from the error dif- 
ference factor, A E f . Most importantly having satellites at high elevation angles with re- 
spect to the termini of the baseline is not detrimental to relative positioning as high 
elevation angles are for point positioning. 

A second area of investigation is the path dependent errors. The ionospheric and 
tropospheric errors, approximated by cosecant functions of the elevation angle, are 
minimized at high elevation angles. Simultaneous high elevation angles for both stations 
is another by-product of high PDOP. 

Finally, the effects of differencing must be explored. Qualitatively, when the satel- 
lites are bunched close together at high PDOP (Figure 3), some of the correlation in 
errors that was lost in lengthening the baseline is restored. The change in tropospheric 
and ionospheric errors with time will also be minimized because: 



d esc e 
dt 



— cos e 

■ 2 
sin e 



d e 
dt 



(6.5) 



and high elevation angles minimize Equation (6.5). Such time derivative minimization 
of errors, will greatly affect the triple difference because the triple difference is itself a 
time difference. 

The sum of the elevation angles for all satellites at both stations (Table 22) reaches 
a maximum after 70 minutes of observations, which is ten minutes later than the PDOP 
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peak. The minimum variance is achieved in the 80 to 90 minute Case 1 time. The run- 
ning mean of the sum of the elevation angle peaks at the 100 minute observing time 
while the running mean elevation angle reaches a local maximum in the 80 to 100 minute 
range. The running means provide a better estimator of the behavior of the accuracy 
and precision because of the accumulation of observations. 



Table 22. ELEVATION ANGLES 



Observation 

Time 


Number 
of SV's 


Sum of 
Elevation 
Angles 


Mean 

Elevation 

Angle 


Running 
Mean of 
Sum 


Running 

Mean 

Elevation 

Angle 


0 


S 


404 


50.5 


404 


50.5 


10 


8 


446 


55. S 


425 


53.1 


20 


10 


509 


50.9 


453 


52.3 


30 


10 


525 


52.5 


471 


52.3 


40 


10 


547 


54.7 


486 


52.8 


50 


10 


559 


55.9 


49S 


53.4 


60 


10 


569 


56.9 


508 


53.9 


70 


10 


562 


56.2 


515 


54.2 


80 


10 


557 


55.7 


520 


54.4 


90 


10 


54S 


54.8 


523 


54.4 


100 


10 


535 


53.5 


524 


54.3 


110 


10 


517 


51.7 


523 


54.1 


120 


8 


505 


50.5 


519 


54.7 


130 


8 


464 


5S.0 


516 


54.9 


140 


8 


452 


56.5 


511 


55.0 


150 


8 


441 


55.1 


507 


55.1 


160 


8 


418 


52.3 


502 


55.0 


170 


8 


396 


59.4 


496 


54.8 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Following a recommendation by Remondi [1984] that more testing of the accuracies 
of triple difference carrier phase measurements was needed, I studied the optimized GPS 
observation times required to achieve geodetic accuracies by partitioning observation 
periods for a 1230-km baseline into ten-minute segments and changing the length and 
starting time of observations. Accuracy was determined by comparing measured 
baseline values with reference values obtained by locating the ends of the baselines from 
very precise VLB I horizontal control points using double difference GPS carrier phase 
measurements over short baselines. 

The 1230-km baseline was directly measured for 28 days over a period of eight weeks 
using Trimble 4000SX receivers. The Trimble supplied Trimvec software was used to 
process the carrier phase measurements. Using broadcast ephemerides in a fixed orbit, 
triple difference carrier phase solution with no ionospheric corrections and a 
tropospheric delay model that used only surface meteorological data, 1.0 part per million 
(ppm) accuracy in the slope distance was achieved for any observing period with a day- 
to-day repeatability better than 1.0 ppm (2a). The AX, AY, and A Z (the components 
of the baseline parallel to the WGSS4 axes) achieved accuracies better than 1O.0 ppm for 
any observing period. AY and AZ repeatabilities were better than 10.0 ppm for any 
observing period, while AX required 30 minutes of observations to reach 10.0 ppm. 

I had not expected to achieve 1.0 ppm accuracy in the slope distance measurements 
because Remondi [1984] had suggested that 1.0 ppm required dual frequency measure- 
ments using highly accurate orbit information and water vapor measurements, and my 
solutions are further burdened by the correlations of the triple differences. The 1.0 ppm 
accuracies were on the order of the uncorrelated single frequency results of previous long 
baseline GPS surveys: Cannon et al. [1985], Bock et al. [1984], Goad et al. [19S5], and 
Mader and Abell [1985]. 

The unexpected 1.0 ppm accuracy can be attributed to low ionospheric activity be- 
cause of the orientation of the baseline, the time of year, and the minimum in the 
sunspot cycle and solar activity during the observing period. Ephemerides errors are 
considered low because of the low age of the ephemerides. Those ameliorating factors 
must be considered before applying the results of these case studies to planning surveys 
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on dillerent baselines at dilTerent times of day and year. Ionospheric errors were on the 
order ol the tropospheric errors. The superposition of the opposing ionospheric (short- 
ening) and tropospheric (lengthening) errors reduced mean slope distance errors. 

Slope distance errors and observation periods were reduced when GPS observations 
started near the infinite, symmetric PDOP peak. The reduction of errors about the 
PDOP peak can be attributed to the simultaneous minimization of: range errors, the 
projection ol range errors onto a baseline, satellite and epoch separation. A high sum 
of the satellite elevation angles for a single station can also be used to determine favor- 
able observing times, as simultaneous high elevation angles correspond to high PDOP 
values. 

My results reconfirm Trimble Navigation's [1987b] recommendation to include the 
infinite PDOP peak in the observing session; however, I could not confirm the Federal 
Geodetic Control Committee's [1986] proposal to stop at a GDOP of 5.0. I could not 
identify any consistent observing stop point that improved accuracy other than to stop 
when four satellites were no longer available. 

When less than 30 minutes of data are used, the goodness of the C A code has a 
great effect upon the carrier phase solution because C/A code positions are used to es- 
timate the carrier phase solution and to compute the time tags of the carrier phase 
measurements. The number of rejected triple difference measurements provided the best 
indicator of the quality of the C/A code estimates and the carrier phase solutions. In- 
cidences of poor accuracy of the carrier phase solution were found to be caused by poor 
C A code estimates. The carrier phase solutions could then be improved by correcting 
the C.'A code solutions, and in turn the carrier phase time tags, by pre-processing the 
C A code measurements or accepting more carrier phase measurements. Pre-processing 
the C A code measurements did not improve the accuracy of the carrier phase solutions 
when the number of rejected triple differences was less than 10% or when the solution 
standard deviation of the slope distance was less than 10 m. 

The long-term average of the best pseudorange solutions produced results compa- 
rable to the mean carrier phase solution, except that the pseudorange solutions had 
variances several times the magnitude of the carrier phase solutions. 

B. RECOMMENDATIONS 

The results of these case studies should be tested on surveys conducted with 
baselines of varying orientation and length and with different satellite configurations. 



51 



Currently, TRIMVEC does not support the use of precise ephemerides or the com- 
putation of uncorrelated triple differences. Should those capabilities be implemented in 
the future, or provided by third-party software, the data should be reprocessed and an- 
alyzed to isolate the tropospheric and ionospheric error contributions and to improve 
the baseline component results. Both of those capabilities would allow better insight 
into the effects of geometry and the tropospheric and ionospheric errors. 

The usual role of the correlated triple difference is to aid in fixing cycle slips when 
the receiver loses lock and to provide the initial estimates for the double difference 
least-squares processing. The current data should be reprocessed in the double differ- 
ence mode to determine whether the correlated triple difference solution provides suffi- 
cient accuracy to fix the cycle slips and estimate the initial integer ambiguities in long 
baseline surveys. 

The standard values used to estimate accuracy could be improved by a more rigor- 
ous fixing of the antenna locations. Both sites should be subjected to a network ad- 
justment, either locally or simultaneously (possibly in conjunction with the NGS VLB! 
crustal motion studies). The extra baselines required for a network solution could be 
measured using the third Naval Postgraduate School receiver. 

Studies should be conducted to determine the effects of using meteorological obser- 
vations far removed from the Sand Point antenna site. A temporary meteorological 
station could be set up near the antenna site (possibly in cooperation with the Weather 
Service Forecast Office (WSFO) located on the grounds of the Western Regional Center 
or with the nearby University of Washington). Future Naval Postgraduate School GPS 
surveys would be aided by portable meteorological instruments, preferably that made 
digital records of the temperature, humidity, and pressure. 

Because of the importance of elevation angles to accurate results, studies should be 
conducted to determine the performance of the TRIMVEC supplied tropospheric mod- 
els( modified Hopfield and Marini ) at various elevation angles. 

While this study has shown the utility of the single-station PDOP, a more complete 
DOP may provide a better satellite selection aid. Such a complete DOP ( or Dilution 
Of Relative Position (DORP)) should incorporate covariances for baseline components, 
satellite orbital errors, receiver timing errors, ionospheric and tropospheric delays, un- 
certainties in the reference station coordinates, and cross-covariances. The DORP 
would also be formulated with respect to the type of differencing to be used. 
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Lastly, the amount of data processed for these ease studies was limited by the speed 
of the microcomputer and the necessity of transferring the results to the Naval Post- 
graduate School mainframe for analysis. As more GPS surveys will be conducted in the 
future, the demands for processing power will increase. It will become imperative that 
the processing software be ported to the mainframe or to a networked mini-computer. 
Those computers will allow multiple access to the software as well as speeding the 
processing. As an interim measure, the RAM of the current ensemble of microcomput- 
ers should be increased to several megabytes. This increase in RAM will allow the use 
of RAM -disks which will speed up the data loading - which is the most time consuming 
of the processing procedures. 
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APPENDIX. BATBLD.BAS LISTING 



10 REM PROGRAM BATBLD.BAS BY R. BOUCHARD NOV87 

20 REM PROGRAM BATBLD: BUILDS THE BATCH FILE FOR TRIMVEC PROCESSING. 

30 REM READS WX OBS FROM SEA. MET FILE 

40 REM DAY 283 REQUIRES ITS OWN BATBLD PROGRAM. 

50 INPUT "batch file name";0FL$ 

60 SP$=" " 

70 Sl$=" command /c tbf nodd. tern " 

80 OPEN OFL$ FOR APPEND AS V 2 

90 DIM A $ ( 18) , JD$( 29 ) ,MM$( 29) , DD$( 29 ) , PH( 29 ) , PM( 29 ) 

95 REM initialize Julian Day array 

100 JD$(n="272": JD$( 2)="273": JU$( 3)="274": JD$(4)="275" 

: JD$(5)-276" 

110 JD$(6)=”279": JD$( 7 )="280": JD$( 8)="28l": JD$(9)="282" 

: JD$(10)="287" 

120 JD$(in="290": JD$( 12)="293": JD$( 13)="294": JD$( 14)="295" 

: JD$( 15 )- 297" 

130 JD$(16)="300": 

140 JD$(20)="308" 

145 JD?(23)="315" 

147 JD$(26)="325": 

150 MM$(1)="09": MM$(2)="09" 

160 FOR 1=3 TO 18: MM$(I)="l0": NEXT I 
170 FOR 1=19 TO 28: MM$(I)="ll": NEXT I 

180 DD$(n="29": DD$(2)="30": DD$(3)="0l": DD$(4)="02": DD$(5)="03" 

: DD$( 6)="06": DD$( 

7)="07": DD$(8)="08": DD$(9)="09": DD$( 10)»"l4": DD$(ll)="l7" 
DD$(12)="20" 



300": 


JD$( 17)=" 301": 


JD$(18)="302": JD$( 19)="307" 


308": 


JD$(21)="309": 


JD$(22)="314" 


315": 


JD$(24)="317": 


JD$(25)="318" 


325": 


JD$(27)="329" 


: JD$(28)="32l" 



DD$(14)="22" 

DD$(19)="03" 

DD$(24)="13" 

DD$(27)="25" 



DD$( 15)="24": DD$( 16)="27": DD$(17)="28" 
DD$(20)="04": DD$( 21)="05": DD$(22)=”lO" 
DD$(25)="14" 

DD$(28)="17" 



A$(2)="a02": A$(3)="a03": A$(4)="a04": A$(5)="a05" 



190 DD$(13)="2l" 

200 DD$(18)="29" 

205 DD$(23)=" ll" 

207 DD$(26)="2l" 

210 A$(n='aOl": 

: A$( 6)="a06" 

220 A$(7)="a07": A$(8)="a08": A$(9)="a09": A$( 10)=”al0": A$(ll)="all" 

: A$(12)="al2" 

230 A$(13)="al3": A$( 14)="al4": A$( 15)="al5": A$(16)="al6" 

: A$( 17)="al7": A$(18)="al8" 

240 FOR 1=1 TO 7: PH(I)=19: NEXT I 
250 FOR 1=8 TO 13: PH(I)=18: NEXT I 

PH( I ) = 1 7 : NEXT I 
PH( I ) = 16 : NEXT I 
PH( I )=15: NEXT I 

280 PM( 1)=34: PM(2)=31: PM(3)=26: PM(4)=23: PM(5)=18: PM(6)=5: PM(7)=1 
290 PM(8)=58: PM(9)=54: PM(10)=33: PM(11)=22: PM(12)=8: PM( 13)=3 
: PM( 14)=59 

300 PM( 15 )=51: PM(16)=39: PM(17)=35: PM(18)=31: PM(19)=14 
310 PM( 20)=6: PM(21)=2: PM(22)=42 

315 PM( 23) =3 7 : PM(24)=29: PM(25)=24: PM(26)=56: PM(27)=39: PM(28)=8 
317 REM begin building the batch file for each day. 

320 FOR IK=1 TO 28 
330 OPEN "i", 1, "sea. met" 



260 FOR 1=14 TO 21 
270 FOR 1=22 TO 25 
275 FOR 1=26 TO 28 
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340 INF$="sa"+JD$(IK)+" ma"+JD$(IK) 

350 S2$=S1$+INF$ 

360 PRINT "FOR JD: ";JD$(IK) 

370 INPUT "start hour"; SH 
380 INPUT "start minute"; SM 
390 PRINT JD$ , PH( IK) , PM( IK) 

400 SiI$=STR$(SH) 

410 SM$=STR$(SM) 

420 PC=SM/60 
430 1=1 

435 REM Find the WX OB 

440 INPUT# 1 ,ID$,IH$,P1(I),T1(I),R1(I),P2(I) ,T2(I),R2(I) 
450 IF( ( ID$=JD$( IK) ) AND( SH=VAL( IH$ ) ) ) THEN 470 
460 GOTO 440 

470 P 1 ( I )=P1( I) -2 : T2(I)=(T2(I)-32)*5/9 

480 1=1+1 

490 INPUT# 1 , ID$,IH$,P1(I),T1(I) ,R1(I) ,P2(I) ,T2(I) ,R2(I) 
500 IF( ID$<>JD$( IK) ) THEN 540 
505 REM Subtract 2 mb from SEA ob 
510 P 1 ( I)=P1( I) -2! 

515 REM Convert NPS Temp to Celsius 
520 T2( I)=(T2( I ) -32)*5/9 
530 GOTO 480 
540 IE=I-1 

550 REM interpolate WX OBS TO GPS START TIME 

560 PR1( 1)= Pl( l)+( Pl( 2) -Pl( 1) )*PC 

570 PR2( 1)=P2( l)+( P2( 2 ) -P2( 1))*PC 

580 TR1( 1 ) =T 1 ( 1)+(T1( 2) -Tl( 1))*PC 

590 TR2( 1)=T2C 1)+(T1( 2) -Tl( 1))*PC 

600 RR1(1)=R1(1)+(R1(2)-R1( 1))*PC 

610 RR2(1)=R2(1)+(R2(2) -R2( 1) )*PC 

615 REM Compute running mean of first hour 

620 PR1( 2)= ( PR1( 1)+P1( 2) )/2 

630 PR2(2)=(PR2( l)+P2(2))/2 

640 TRl(2)=(TRl(l)+Tl(2))/2 

650 TR2( 2)=(TR2( 1)+T2( 2) )/2 

660 RRl(2)=(RRl(l)+Rl(2))/2 

670 RR2( 2)=(RR2( 1)+R2( 2 ) )/2 

680 WT=1-PC 



685 REM Compute remaining running means 



690 FOR 1= 3 TO IE 
700 IA=I - 1 
710 IWT=1+WT 
720 PR1( I)=( PR1( IA)*WT 
730 PR2(I)=(PR2(IA)*WT 
740 TR1( I)=(TR1( IA)*WT 
750 TR2( I )=(TR2( IA)*WT 
760 RR1( I)=(RR1( IA)*WT 
770 RR2(I)=(RR2(IA)*WT 
780 WT=IWT 
790 NEXT I 



+ ((Pl(IA)+Pl(I))/2))/( IWT) 
+ ( ( P2( IA)+P2( I))/2))/( IWT) 
+ ((Tl(IA)+Tl(I))/2))/( IWT) 
+ ( (Tl( IA)+T2( I))/2))/( IWT) 
+ ((Rl(IA)+Rl(I))/2))/( IWT) 
+ ( (R2( IA)+R2( I))/2))/( IWT) 



800 M1=SM 



810 H=SH 
820 C=0 
830 1=2 

840 DHR=PH( IK) -SH: DMIN=PM( IK) -SM 
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850 IF(DMIN<0) THEN DHR=DHR- 1: DMIN=DMIN+60 
860 NP=DHR*6+( DMIN/ 10 ) 

865 REM Interpolate hourly running means to 10 minute intervals 
870 FOR IL=1 TO NP 
880 M1=M1+10 

890 IF(M1>60) THEN Ml=Ml-60: H=H+1: 1=1+1 
900 WT=Ml/60 

910 IF(M1=60) THEN Ml=Ml-60: H=H+1: 1=1+1 

920 MP=M1 

930 HP=H 

940 IWT=1-WT 

950 IA=I - 1 

960 PW1=PR1(IA)*(IWT)+PR1(I)*WT 
9 70 PW1=INT( ( 10*PW1)+. 5)/ 10 
980 PW2=PR2( IA)*( IWT)+PR2( I )*WT 
990 PW2=INT((10*PW2)+. 5) / 10 
1000 TW1=TR1( IA)*( IWT)+TR1( I )*WT 
1010 TW1=INT(( 10*TW1)+. 5)/10 
1020 TW2=TR2( IA)*( IWT)+TR2( I )*WT 
1030 TW2=INT((10*TW2)+. 5)/10 
1040 RW 1=RR 1 ( IA)*( IWT)+RR1( I )*WT 
1050 RW1=INT((10*RW1)+. 5)/10 
1060 RW2=RR2( IA)*( IWT)+RR2( I )*WT 
1070 RW2=INT( ( 10*RW2)+. 5)/ 10 
1080 REM convert to strings 
1090 PS1$=LEFT$(STR$(PW1) ,7) 

1100 PS2$=LEFT$( STR$( PW2) ,7) 

1110 TS1$=LEFT$(STR$(TW1) ,5) 

1120 TS2$=LEFT$(STR$(TW2) ,5) 

1130 RS1$=LEFT$( STR$(RW1) ,5) 

1140 RS2$=LEFT$(STR$(RW2) ,5) 

1150 HS$=STR$(HP) 

1160 S3$=S2$+SP$+A$( IL)+". ”+JD$(IK) 

1170 MS$=STR$(MP) 

1180 PRINT S3$+PS1$+TS1$+RS1$+PS2$+TS2$+RS2$+SP$+MM$(IK)+SP$ 
+DD$( IK)+HS$+MS$+SH$+SM$ 

1190 PRINT//2, S3$+PS1$+TS1$+RS1$+PS2$+TS2$+RS2$+SP$+MM$(IK) 
+SP$+DD$( IK)+HS$+MS$+SH$+SM$ 

1200 LPRINT S3$+PSl$+TSl$+RSl$+PS2$+TS2$+RS2$- i -SP$+MM$( IK)+SP$ 
+DD$( IK)+HS$+MS$+SH$+SM$ 

1210 NEXT IL 
1220 CLOSE , U 
1230 NEXT IK 
1240 CLOSE #2 
1250 END 
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